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AESTRACT 


The  objectives  of  this  paper  are  two-fold.   In  theory  developed 
previously  it  is  shown  how  to  express  any      y  state  electromagnetic  field 
in  the  form  of  an  asymptotic  series.   The  "basic  variable  L  -  this  series  is 
the  wave  length,  i.e.,  the  first  term  is  given  for  zero  .  ;.-.   .ength.  The 
present  paper  serves  first  as  an  illusv     n  of  how  th        on  asymp- 
totic expansions  can  "be  applied  to  physic     Dhlems.  Sec       the  appli- 
cation itself  is  of  interest  and  importance  i:      -enerai        f  diffrac- 
tion through  lenses  and  as  such  it  discusses  a  pre     now  ac     ■   pursued 
in  current  research.   Insofar  as  application  to  diffraction  is  ec  earned 
this  paper  gives  a  wave  formulation  of  ,       at  generalized  pre ._em  . f 
diffraction  through  lenses  and  obtains  e  solution   f  the  for     .en  in 
terms  of  integrals.  At  this  point  the  theory  of  asymptotie       .ons      s 

applied  to  evaluate  the  integrals.  'It   turns  out  that  the  leading  tcr 

the  asymptotic  expansions  of  the  exectrci — gnetic  fields,  i.e.,  the  ter 
which  represent  the  "behavior  of  the  field  either  for  very  small  wave  len^  - 
at  finite  distances  from  the  lens  or  for  arbitrary  wave  length  at  infinite 
distances  from  the  lens  are  precisely  the  terms  given  by  geometrical  optics. 
The  higher  order  terms  of  the  expansions,  though  not  evaluated  explicitly 
in  this  paper,  give  the  diffraction  effects  which  wave  theory,  as  opposed  to 
geometrical  optics,  encompasses. 
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1.  Introduction 

1.1.  In  an  earlier  paper  the  author  showed  how  the  complex  spatial 
factors  u  and  V  of  any  steady  state  electromagnetic  field  can  each  be  repre- 
sented "by  an  asymptotic  series  in  which  the  "basic  variable  is  the  wave-length 
The  coefficients  of  the  terras  in  each  asymptotic  representation  depend  primar- 
ily upon  the  discontinuities  with  respect  to  time  in  what  is  called  the  pulse 
solution  of  Maxwell's  equations.  This  solution  is  the  field  set  up  by  charges 
which  have  the  same  geometrical  distribution  as  the  source  of  the  steady  state 
field  but  which  are  suddenly  allowed  to  act. 

The  present  paper  shows  how  this  general  theory  can  be  applied 
to  the  evaluation  of  integrals  (see  formulas  (U.15)  v/hich  arise  in  the  study 
of  diffraction  of  electromagnetic  waves  through  a  medium  of  variable  electro- 
magnetic parameters  £,  JUL,   and^".   The  mathematical  formulation  of  the 
diffraction  problem  herein  treated  is  a  generalization  of  the  common  optical 
problem  of  diffraction  through  a  lens.  The  formulation  of  the  solution  to  the 
diffraction  problem  is  not  in  itself  especially  novel.  It  leads  in  fact  to  a 
well  known  type  of  diffraction  integral  which  is  unfortunately  too  complicated 
to  be  readily  evaluated  except  in  the  simplest  of  cases.  The  essence  of  this 
paper  will  then  consist  in  showing  hovr  '. tegrals  of  this  type  can  be  approxi- 
mately evaluated  by  the  use  of  the  general  theory  of  asymptotic  expansions 
referred  to  above.   It  will  be  seen  too  that  this  method  of  evaluation  is 
closely  related  to  the  stationary  phase  method  as  applied  to  integrals. 

Unfortunately  even  the  diffraction  integrals  herein  presented 
as  the  solution  to  the  problem  of  the  diffraction  caused  by  a  medium  of  vari- 
able electromagnetic  parameters  are  not  exact  solutions.  Rather  the  diffracted 
field  sought  is  required  to  satisfy  some  reasonable  physical  conditions  and  the 
integrals  are  shown  to  meet  these  conditions. 

While  one  purpose  of  this  paper  is  to  illustrate  the  application 
of  the  general  theory  of  asymptotic  expansions,  another  is  the  development  of 
the  diffraction  theory  of  optical  aberrations.  The  physical  meaning  and 


1.  Luneberg,  R.X. :  Asymptotic  Expansion  of  Steady  State  Electromagnetic  Fields, 
Research  Report  Ho.  3M-1U,  submitted  by  Hew  York  University  to  the  Geo- 
physical Research  Directorate  of  the  Air  Force  Cambridge  Research  Labora- 
tories, AMC,  U.S. Air  Force.  A  more  detailed  review  of  this  earlier  theory 
will  be  given  in  article  2.  Report  No.  EM-lU  will  be  referred  to  as  the 
first  paper. 
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p 
importance  of  such,  theory  is  discussed  "by  Zerrnike. 

An  analogous  asymptotic  evaluation,  of  diffraction  integrals  is 

obtained  by  Van  Kampen. ^    3y  using  Kirchoff 's  principle  Van  Kampen  expresses 

an;/  component  of  the  electric  or  magnetic   field  in  the  form  of 

Jj-q     kg(x,y;e       ^X*      dxdy    where,    in  practical  cases,   D  can  he  the  domain 

of  the  aperture  in  a  "black  screen  provided  D  is  hounded  by  one  or  more 
analytic  curves.  The  asymptotic  behavior  of  this  integral  for  large  k 
depends,   Van  Kampen  finds,   on  the  behavior  of  the  integrand  at   three  kinds 

of  ■ooints;      those  at  which  the  exponent  is  stationary,    i.e.;     l-— i-  =  -£—    =  0: 

'       gx   ®y 

those  at  which  the  exponent  is  stationary  with  respect  to  values  on  the 

boundary  of  D;   and  corner  points,  or  points  where  the  boundary  of  D  fails 

to  be  analytic. 

These  same  types  of  points  appear  in  the  theory  of  this  paper, 
though  the  method  of  obtaining  the  asymptotic  expansion  is  different.  Van 
Kampen  and  the  present  author  find  that  points  of  the  second  and  third  kind 
do  not  contribute  to  the  zero  order  term  in  k.  Hence  their  contribution 
represents  diffract ive  effects  as  opposed  to  geometrical  optics. 

The  diffraction  study  of  this  paper  is  closely  related  to 
numerous  other  studies  of  diffraction  through  apertures.  For  apertures  in 
infinitely  thin  screens  quite  different  methods  have  been  employed  and, 
even  when  approximate,  would  seem  to  be  more  useful.  For  thick  screens, 
however,  the  method  of  this  paper  provides  an  approach  which  can  be  em- 
ployed if  the  geometrical  optics  solutions  is  known. 

It  should  be  noted  that  while  the  general  .~~j?y   of  asymptotic 
development  of  steady  state  fields,  as  presented  in  the  paper  of  footnote  1, 
is  not  restricted  to  fields  in  homogeneous  media  the  application  made  to  the 
diffraction  theory  in  this  paper  is  so  restricted.  The  reason  for  this  is 
that  the  integral  formulation  of  the  diffraction  problem  is  obtained  by  the 
use  of  Green's  theorem  and  since  use  of  this  theorem  calls  for  knowledge  of 
a  suitable  Green's  function  for  the  space  of  the  diffracted  field,  present 
knowledge  restricts  the  use  of  the  theorem  to  homogeneous  media. 

2.  Zerrnike,  F. :  Diffraction  and  Optical  Image  Formation,  Physical  Soc. ,?roc. , 
v.  61,  1948,  pp.  158-164. 

3.  Van  Kampen,  N.G.:  An  Asymptotic  Treatment  of  Diffraction  Problems,  Physica, 
V.  14,  Jan.  1949,  PP.  575-5S9. 

4.  See,  for  example,  Levins,  Harold  and  Schwinger,  Julian:   On  the  Theory  of 
Diffraction  by  an  Aperture  in  an  Infinite  Plane  Screen.I.  Physical  Review, 
V.74,  Oct. 15,  194S,pp.95S-974.  The  Introduction  to  this  paper  describes 
other  methods  of  attack  on  the  problem  of  diffraction  through  an  aperture 
in  an  infinitely  thin  screen. 
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It  is  hoped  that  the  general  theory  of  asymptotic  expansion 

of  steady  state  fields  will  find  application  to  other  integrals  arising  in 

physically  quite  different  problems  which  nevertheless  lead  to  integrals 

such  as  arise  in  the  diffraction  problem.  One  such  problem  is  that  of 

electromagnetic  propagation  in  non-homogeneous  media.  At  the  present  time, 

however,  successful  application  of  the  asymptotic  theory  to  the  integrals  of 

the  propagation  problems  has  not  been  achieved. 

1.2.  Because  the  material  of  this  paper  is  both  lengthy  and  detailed 
a  summary  will  be  presented  in  this  introduction. 

The  x,y,z  space  is  divided  into  three  domains,  the  domain 
z  <  z  ,  called  the  object  space,  the  domain  z  >  z.  called  the  image  space 
and  the  domain  z  <  z  <  z1  in  which  the  material  constants  <£  ,  ^4 ,  0~  vary 

with  x,y,z.  We  shall  refer  to  this  domain  as  the  lens  and  to  the  entire 

object  and  image  space  and  lens  as  the  instrument.  Vfe  assume  that  £.  =   n  in 

2  2      2 

the  object  and  £  =  n  in  the  image  space,  n  and  n  constant;  we  assume  also 

that/^  =  1  and  (J~-   0  in  object  and  image  space.  A  time  periodic  point  source 

is  located  at  a  point  P  of  the  object 


space.  She  problem  is  to  find  the  time 
periodic  electromagnetic  field. 


E=  Real  part  of  "u(x,y,z)e      and 

E  =  Real  part  of  v(x,y,z)e~XCOt   which 
-.     ',  is  established,  in  the  image  space, 

under  the  influence  of  this  source.  The  vectors  "u  and  V*  are  complex  functions 
a(x,y,z)  +  ib(x,y,z)  and  c(x,y,z)  +  id(x,y,z)  respectively.  The  vectors  u"  and 
v* then  are  solutions  of  the  time  free  Maxwell's  equations 


(1.1) 

_  co    _  2rr 


curl  v  +  i  €  ku"  =  0 


curl  u"  -  i  M  lev*  =  0 


/-• 


uhere  k 


c     X 
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The  solution  u,v  of  our  problem  must  "be  related  to  the  geo- 
metical  optics  of  our  instrument,  in  particular  to  the  system  of  light  rays 
which  originates  at  the  point  source  P  .  This  relation  must  he  such  that 

for  small  wavelengths ,X,  that  is,  for  large  values  of  k  =  i=H  ,    the  vectors 
U,v  show  an  asymptotic  "behavior  which  expresses  the  geometric  optical  situa- 
tion created  "by  the  instrument.   In  other  words  we  are  concerned  with  a  cer- 
tain prohlem  of  correspondence:  To  find  a  solution  u,v"  of  (l.l)  which,  for 
k  -»  oo,  shows  a  prescribed  asymptotic  "behavior. 

We  characterize  the  rays  of  the  instrument  by  Hamilton's  mixed 
characteristic  W(xo,yo>z  ;  ptq).  This  is  a  function  of  the  coordinates 

x  ,y  , z  of  the  source  and  of  the  direction  numbers  p,q,r  of  the  rays  in  the 

image  space.  It  determines  the  optical  length  of  the  ray  between  the  point  P 

o 

and  the  base  point  Q,     of  the  perpendicular  dropped  from 

the  origin  0  onto  the  ray  p,q,r.     T,.'a  can  determine  all  v G>tc;>r)       q, 

essential  characteristics  of  the  ray  system  with  the  p 

aid  of  tf. 


c*ty,z.) 


Figure. 2 
If  x,y,  z  is  any  point    in  the  image  space  and  Q,  the  base  of 
the  perpendicular  from  x,y,z  to  the   ray  with  direction  numbera  p,q,r  then 
the  function 

(1.2)  y~,y,z)  =  px  +  qy  +  rz  v  './(p,q) 

1 

represents  the  optical  length  of  the  path  from  P  to  Q,,  that  is,  the  optical 
length  P  Qo  plus  the  optical  length  &  Q.   If  now  we  determine  the  ray  p,q,r 

which  passes  through  a  given  point  x,y,z,  express  p,q,  and  r  in  terms  of  x,y,z, 
and  substitute  these  values  of  p,q,  and  r  in  0  we  obtain  the  function  >^"(x,y,z) 
which  gives  the  optical  length  from  P  to  (x,y,z ).  The  surfaces "^(x.y,  z)  = 

constant  are  the  wave  fronts  of  the  system,  that  is,  the  surfaces  orthogonal  to 
the  rays. 

There  may,  of  course,  exist  several  wave  fronts  through  the  same 
point  x,y,z,  Just  as  there  may  exist  several  rays  which  connect  this  point  with 
the  source  P 
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Let  &CO  be  the  surface  element  in  which  a  narrow  "tube"  of 
rays  intersects  a  wave  front  >^"  =  constant.  We  measure  the  expansion  of  the 
ray  system  by  the  expression 


do)V  1  -  p  -  0/ 

which  in  our  case  of  a  homogeneous  space  is  equal  to  the  Gaussian  curvature 
of  the  wave  front.  If  the  light  distribution  along  a  ray  is  observed  then, 
according  to  geometrical  optics,  the  polarization  of  light  is  unchanged 
whereas  intensity  and  energy  flux  vary  proportional  to  the  above  quantity  K. 
We  may  describe  these  observations  mathematically  by  a  vector  f ieldtj(x,y,z); 
the  direction  of  these  vectors  determines  the  polarization  of  the  light,  the 
absolute  value  |U|  ,  the  energy  and  flux.  The  vectors  U  along  a  fixed  ray 
p,q,r  are  normal  to  the  ray  and  are  related  "by  an  equation  tf  =  X(p,q)"v...' 
where  the  vector  A(p,q;  is  constant  alor.-;  the  ray. 

The  actual  solution  of  the  diffraction  problem  requires  the 
integration  of  Maxwell's  equations  (l.l)  where  in  the  middle  medium  £  and/U 
are  given  functions  of  x,y,z.   On  account  of  the  math^_.;_cal  difficulties 
which  oppose  such  an  approach  one  prefers  to  establish  solutions  in  the 
image  space  by  making  reasonable  hypotheses: 

1)  The  vectors  u  and  v  are  analytical  solutions  of  (l.l)  which  are 
regular  and  uniformly  bounded  in  the  to-:al  xsy,z  space. 

\  u    "»ik  *\^*"  —£> 

2)  The  vectors  — — -  e      attain  definite  limit  values  A(p,q) 

and  B(p,q;  respectively,  if  infinity  is  approached  along  a  ray  p,q,r  of  the 
instrument.  In  other  words  it  is  assumed  that  at  infinity  both  the  actual 
electromagnetic  field  and  the  geometrical  optics  approximation  have  the  same 
asymptotic  behavior. 

By  application  of  Green's  theorem  it  is  shown  that  a  solution  ■ 
satisfying  these  hypotheses  can  be  obtained  by  superposition  of  plane  waves, 
namely  by  the  integrals 


ff  A(p,q)    eik0 


->      ik 

u  =  -  —  JJ  "^'^-  e"  '  dpdq 

C         D  •*/!   2 2 

Vn  -  p  -  q 


(1.3) 


k-/f  ~5(p,e)    Qik0 


D  -t/~2   2  2 

v  n  -p  -  q 
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The  domain  of  integration,  D,  is  the  range  of  all  p  and  q  values  of  rays  which 

cover  the  image  space  and  which  are  therefore  determined  by  the  medium  "between 

zq  and  z1  or  the  lens.   It  is  a  domain  of  the  p,q  plane  which  lies  inside  the 

2    2    2 
unit  circle  p  +  q  =  n  .  The  integrals  (1.3)  are  solutions  of  the  equations 

(l.l)  for  any  choice  of  vectors  A,  B  provided  that  the  conditions 

p  x  B  <■  A  =  0 
px;^  -  B  =  0  , 

wherein  p*  is  the  vector  (p,q,r),  are  satisfied. 

Evaluation  of  the  integrals  (1.3)  by  known  functions  is  possible 
only  in  a  few  simple  cases.  Thus  we  are  led  to  the  problem  of  determining  the 

principal  terms  of  these  integrals  by  an  asymptotic  development  with  respect 

2n 

to  a  parameter  proportional  to  the  quantity  k  =     —  .     The  solution  of  this 

A 

problem  then  should  provide  the  answer  to  the  question  whether  the  integrals 
solve  the  above  formulated  problem  of  correspondence  to  a  prescribed  situation 
of  geometrical  optics.   The  answer  will  be  in  the  affirmative:  'The  -orinci-oal 
terms  of  the  integrals  (I.3)  are  given  by  the  expressions  u"  =  A  "—-  e    , 


v  =  3  __ 
o      k 


_1:  ■> 


Me   approach  the  problem  of  asymptotic  development  as  follows: 
We  introduce  the  vectors 


P(x,y,z,t)  = 


2TTC 


(l.U) 


Q,(x,y,z,t;  = 


2TTC 


Figure  3 


(1.5) 


u 


// 

(b  <co 

fl 


),q) 


-»/2        2     2 

v  n  -  p  -q 

3*(p,q) 

-i/~2        2     2 

V n  -  p  -q 


dpdq 


dpdq 


The  domain  of  integration  is  that  part  of 
the  domain  D  where 


4>=   px  *   qy  +  rz  +  V  <  ct  . 


'tie  then  may  write  the  multiple   integrals   in 

the  form  of  simple  integrals: 
•  <=» 


=  -ico  J     P»(t 


>,   icut      .. 

je  dt 


0 
00 


v  : 


-io,  /'  ti(t)eia?t     dt 
o 


-s- 
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The  vectors  P'  and  Q,'  are,  at  any  fixed  point  x,y,z,  different  from  zero  only 
in  a  finite  interval  T  <  t  <  T?,  determined  "by  the  extreme  values  of  the  func- 

1  yA  -"      -•■ 

tion  —  p  in  the  domain  D.  In  general  ?  and  Q,  will  depend  analytically  on  t; 
c 

->         -> 

there  are  however  certain  points  t  u  inside  the  ahove  interval  where  P  and  Q, 

or  their  derivatives  are  discontinuous.  These  points  can  he  found  hy  study- 
ing the  topology  of  the  level  lines  (y  -   const  in  D  . 

Consider  for  example  a  point  x,y,z  in  the  regular  part  of  the 
field  of  rays  (that  is,  only  one  ray  passes  through  x,y,z).   In  this  case  cp 
possesses  exactly  one  stationary  point,  P  ,  in  the  interior  of  D;   the  associ- 
ated value,  ct  ,  is  either  a  maximum  or  minimum  of  <p  and  equals  the  optical 

■}£"(x,y,z)  from  the  source  to  the  point  x8y,z.   The  level  lines  £$  =  constant 


thus  surround  the  stationary  point  P  in 


the  manner  indicated  in  Figure  U.  There 
exist  two  level  lines 

^=  ctx  ,    s  =  ct2 

(perhaps  more)  which  are  tangential  to   the 


"boundary  of  D  at  points  P.    and  P_  . 


Ve 


can  _'— „  chase  level  lines  hy  determining 
the  stationary  points  of  the  boundary  values 

of  the  function  (p   .   If  X  and  IT  are  analytic  in  D  then  P'(t)  and  ^"'(t)  will 
he  analytic  functions  of  t  except  at  the  points 


t  = 
o 


"& 


1    c   1 


c  rZ 


given  hy  the  points  P  ,P.  ,P_  of  3?igure  \.     The  components  of  the  vectors 
P* and  "q  thus  are  functions  of  t  of  the  type  shown  in  Figure  5» 


?G> 


■*.£ 


Figure  5 
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It  is  perhaps  clear  from  the  fact  that  P  and  Q  will  experience  sudden  changes 
in  the  rate  of  increase  with  t  as  t  goes  through  the  values  t  ,tn ,  and  t 
that  these  values  of  t  will  "be  points  of  discontinuity  for  ?'(t)  and  Q'(t). 
If  x,y,z  lies  outside  the  field  of  the  rays  then  there  is  no 
interior  stationary  pc   t  of  0  "but  again  two  or  more  stationary  boundary 
points  P  and  P  .  Hence  there  are  only  two  joints 

->     •*• 

where  P1  and  Q'  are  irregular. 

Figure  6 

If,  finally,  x,y,z  lies  in  the  irregular  part  of  the  ray  system, 
that  is,  if  more  than  one  ray  passes  through  x,y,z,  then  there  will  he  more 
than  one  stationary  point  in  the  interior  and  a  corresponding  number  of 
stationary  boundary  points.  In  this  case  there  exist  a  certain  finite  number 

of  points  t^  =  r  N^  (x,y,z)  where  P'  and  Q,'  are  non-analytic. 

Let  us  now  assume  that  P  and  Q,  are  sectionally  analytic  functions 
of  t,  the  different  branches  being  separated  by  a  discrete,  finite  set  of 
points 

**  =  Iv^y.z)   ,  <=  i u. 

The  points  t^  are  points  of  discontinuity  of  P*'(t)  and^'(t). 

The  representation  (1.5)  for  U*  and  v"  permits  us  to  identify  P*' (t)  andcfl(t) 

with  the  pulse  solution  E  (t)  and  H  (t)  mentioned  earlier  and  described  in 

00 

article  2.      It  is  therefore  possible   to  apply  the  results  of  the  asymptotic 

theory  presented  in  the  earlier  report.      It  follows  from  this  theory  that  u" 

and  v  admit  asymptotic  developments   of  the  form 


U     =  <£-T 


^=1     H*  (x,y,z)e 


(1.6) 

N 


▼  =   ^     V  (x,y,z)e      '- 


wherein  the  values  of  the  vectors  U~  and  Y^     depend  upon  the  behavior  of  3  (t) 
and  H  (t)  in  the  immediate  neighborhood  of  t  =  tcc=  —  ">,/&;('x ,V *  7-)« 
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->  .  . 

If,  for  example,  E (tj  allows,  in  the  neighborhood  of  t  =  t^  , 

representation  of  the  following  fora 

CO 

f0u)  =  (t«  -  tr  «—  c^t^  -t)n  ,  t  <  ^ 

(1.7)  <~ 

B0(t)  «  (t  -  toe  )m  ^o  Da(t  -  t«)n  .  t  >  toe. 

then  the  function  XJ-cay  be  developed  into  an  asymptotic  series  of  the  fora 

(1.8)  _   ,   c   Y    K^n*!^)     Gar(n^j)  | 


where  JL  (%)  is  the  Gaama  function. 

The  problem  of  asymptotic  evaluation  of  the  diffraction  integrals 
is  thus  reduced  to  the  problem  of  finding  the  character  of  the  irregularities 
of  the  vectors  E  and  H  at  certain  singular  points  tec  .  These  singular 

points  are  determined  "by  the  system  of  the  level  lines  4p  =  constant  in  D,  in 
particular  by  the  stationary  points  of  the  interior  or  of  the  boundary  values 
of  the  phase  function  <p  %     We  mention  that  an  interior  regular  maximum  or  mini- 
mum point  of  (p  introduces  an  ordinary  discontinuity  (jZ  =  a  =  0)  of  E  (t);   a 
stationary  boundary  point,  however,  introduces  a  square  root  discontinuity 
(  -c  =  1/2  or  a  =  1/2).  This  shows  that  the  asymptotic  contribution  associated 
with  stationary  boundary  points  (?..  and  P_  in  Fig.  h)  are  of  higher  order  in 

r;  than  those  from  interior  stationary  points.  One  also  finds  that  the  first 

coefficient  C  or  D  in  case  J(=>  m  =  0  has  the  value  AT/K  so  that  the  principal 

terms  in  the  development  (1.6)  have  the  form 

(1.9)  t  -  pH   ."*  *o  -  |vf  a1^ 

as  was  stated  above. 

If  the  function  r)  possesses  more  than  one  interior  stationary 
point  then  a  finite  sum  of  terms  (1.9)  is  obtained  as  principal  asymptotic 
contribution.  This  is  the  case  if  more  than  one  ray  passes  through  the  point 
x,y,s  in  question. 
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2.  Hsvic;;  of  General  Theory  on  Asymptotic  Expansion  of  Steady  State  Fields. 

This  article  will  present  without  proofs   the  essential  material 
of  the  theory  of  asymptotic  expansions  previously  developed^     and  applied 
later  in  this  paper. 

-.    the  pulse  solution  or  pulse  field  satisfying  Maxwell's  equa- 
tion we  shall  mean  the  functions  E  (x,y8z,t)  and  H  (s,y,z,t)  whichare  set  up 

"by  a  distribution  of  charge  which  at  time  t  =  0  is   suddenly  allowed  to  act. 
The  charge   itself  is  characterized  hy  a  vector  function  5*  =  g"(x8y  z)  f(t) 
where  f(t)  is  0  for  t   <  0  and  =  1  for  t  >  0.     The  vector  function  3"  character- 
izes the  charge  in  the  sense  that  the  current  set  up  "by  enforced  motion  of  the 

=     -^j    Et(x8y,z.t)  and  the  charge  distribution  itself 

is  given  "by  />  =  -  div  g.     In  the  case  f(t)  is  the  function  just  described,   the 
effect  of  the  sudden  creation  at   t  =  0  of  the  charge  is  to  set  up  an  electro- 
magnetic field  which; as   t  -5>oOj  approaches  the   stationary  electric  field  of  a 
charge  distribution  />  =  -  div  g*.     That  is  E  (x,y,z,t)  - >0  as   t  -. >oo    for  all 
x,y,z. 

If  the  charge  is  characterized  "by  a  function  3?  =  g"(x,y,z)e 
so  that   the  enforced  current  is  periodic  while  the  charge  distribution  is 
/0=  -div  g"then  it  i3  shown  in  the  earlier  paper  that  the  field  set  up  hy 

the  harmonic  charge  oscillations  is,    at  t  =    co   t   E  =  Real  part  of  u  e  and 

-jv  «^    —icn»t 

H  =  Seal  part  of  v  e  where 

•s.  rOO    ».  ^ 

u(x,y8z)  =  So(x,y,z,co)  -  i&>J        [Eo(x,y,z,<r )  -  EQ(xsy, z,  oo)J   e       rdf 

(2.1) 

v(x,y8z)  =  fu.y.z.oo)  -  itof*      [f  (xsy8z,  *0  -  E^x.y,  z,  oo)]    e^^df 

"'o 

wherein  the  term  Eo(x,y,z,oo)  =  Eq(oo)  has  "been  retained  for  symmetry.  The 

relation  (2.1)  shows  that  harmonic  oscillations  of  charge  produce  at  t  =  oo 
a  harmonic  steady  state  field. 

Prom  the  purely  mathematical  relationship  of  "u  and  v"  to  the 

pot 


fields  EQ(x,y,z,t)   and  Ho(x,y,z,t)  as  given  hy   (2.1)  it  is  possible  to  derive 


asymptotic  formulas  for  u*  and  v*  in  terms  of  properties  of  the  functions  if   and 
2     and  the  circular  frequency  Co  of  the  steady  state  field.     One  introduces  the 

5.     See  footnote  1. 
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vector  fields  A  and  B  defined  by 


~fu)  = 


2ni  J 


;r> 


(r-  t)1 


ar- 


(2.2) 


*>  "  *'/" 


so(r) 


dr 


(r-  t)' 


wherein,  the  x8y,z  dependence  of  A,B,E  ,  and  E  are  -understood  and  of  course 
are  fixed.  Also,  t  is  the  complex  variable  t.,  +  it   .  It  is  then  shown 

that  the  singularities  of  A  and  £  as  functions  of  t  are  the  singularities  of 
t     and  H*  ,  respectively,  as  functions  of  the  real  variable  T  .  The  singu- 

larities  of  E  and  E  occur  at  values  ^;  =  —  V^Cx.y.z),  that  is,  they 

0  0  c 

occur  at  values  of  T  which  lie  on  hypersurfaces  1^(x,y,z)  and  may  arise 
physically  from  the  sudden  passage  of  wave-fronts  across   (x,y,z).      These  fronts 
themselves  can  arise  from  either  the  directly  transmitted  wave-front  created 
by  the  sudden  action  of  the  charge  as  the  front  expands  through  space  or  from 
the  passage  of  fronts  reflected  by  or  refracted  through  discontinuities  in  the 
whole  space. 

The  basic  result  on  the  asymptotic  form  of  u  and  V  can  now  be 
At  a  point  x,y,z  for  which  g(x,y8z)  =  0, 


stated. 


(2.3) 


where 


-*-  -■?     ik  y& 

u  =      <£-=■     IL/  e      roc 

cc=l 


*-      T*    t,eik^ 


%.    J  l(2ht   +r)  e1 


car 


dr 


iM 


M 


% 


v  T)   e  d7^ 


The  path  M  surrounds  the  positive  imaginary  axis 
of  the  coa-olex  7-plane.     Again  the  dependence  of 


<r2 


..: 


7 


->     A. 


-> 


ut*    U<     .Vet     ,  A,    and  B  on  x,y,z   is  understood. 


Figure  7 
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Since  the  values  of  U^  and  Y&    depend  upon  the  singularities 

->    ^r  -5-     -> 

of  the  integrands  A  and  B  and  since  the  singularities  of  A  and  B  are  those 

of  E  and  H  with  respect  to  tine,  it.  is  possible  to  determine  a  general 

->>       -N 
result  for  the  Uc4  and  lac    "by  considering  a  general  form  for  the  singularity 

in  E  and  A     .     The  general  result  is  as  follows:  Let  t  =  t^  (=  —  >& )  be 

0       0  c 

a  discontinuity  of  E  (t)and  let  E   have  the  form 

X—        C 
->      i'P(toc-t)  4r0  — r  (t^-t)      ,  for  t  <  t* 
(2.h)    EQ(t)  =) 


St 


-t*)  ^-0  -jy  (t-t<^r      .  for  t  >  t^ 


In  this  expression  C     and  D     are  constants.  P  and  Q,  are  analytic  functions   of 
r  n  n 

t^  -t  which,  however,    are  singular  at  t  =  %ot    .      For  exsziple  ?  and  Q,  may  he 

of  the  form  (t  -  tocT~     where  k  is  fractional  or  may  "be  of  the   form  log(t-t£^  ). 

->• 
If  E  (t)  has  the  form  given  by  (2.4)  in  the  neighborhood  of  a 

0  vkt  ■*•  f-Y,  \  xC£/r 

singularity  t  =  t^  ,    tu  =     ~— '  ■>   then  U<^    „  which  ha3  the  value  _/ A(tc^+7";e  d.'r 

c 

is  given  by 

C2.5)        ~    (t  roo     .  t  r*» 

■f  >        j     n  1  /        -•"/   t     wn  -t    ..  n  1  -,t  s.n  -t 

The  same  result  applied  to  H     and  T$x     • 

-> 

When  one  obtains  the  various  t^t  and-  1L#  for  all  singularities 

^   s  _i_£i    one  has  but  to  substitute  in  the  summation  formulas  of  (2.3) 
to  obtain  the  asymptotic  expressions  for  u  and  v.  The  dependence  of  u  and 
"v"  on  co  or  on  X  results  through  the  expressions  for  TJc{   and  Ycl    &s  may  be  seen 
in  (2.5). 

For  many  forms  of  the  functions  P  and  Q,  the  integrals  in  Ucc  are 

is  of  the  form  t  and  Q,  of  the  form  t 
where  ^  and  m  are  greater  than  -1,  then  the  integrals  are  readily  evaluated 
through  the  Gamma  function.  In  fact  for  these  P  and  Q,, 


i  D  Pdx+Hm)      C  J(n+l+£) 


U*  =  (-W*»     "  (W* 


The  singularity  in  ZQ   at  t-t^  may  be  a  finite  or  infinite  dis- 

continuity,  or  B  may  "be  continuous  with,  a  discontinuity  in  one  of  its 

derivatives.  It  is  apparent  from  (2.4)  and  (2.6)  that  in  the  expansion  of 

Uathe  first  term  (that  is  for  n=0)  withc^-xo  is  the  discontinuity  in 

o 
at  t-ta  .   (If  Eq  is  continuous  with  some  discontinuous  derivative,  the 

discontinuity  in  %Q   is  obviously  0).  Hence  for  large  CJt   the  sum  of  all 

such  first  terms  would  give  the  first  terms  in  the  expansions  of  u"and"v, 

by  virtue  of  formula  (2.3),  The  latter  terms  would  be5~rf*  1  eik^-  ^^ 

<8»"i   ik-  ,  -> 

^al     e     »  where^j  .aenotes  the  values  of  the  discontinuity  in  E  at 

t=t^>  .  The  recult  apparent  here  from  the  somewhat  specialized  result 
holds  generally  . 

3.  The  Diffraction  Problem  of  Optical  In3trun3nts. 

3.1. We  formulate  first  the  general  diffraction  problem  which  this 
paper  treats.  The  problem  is  a  generalisation  of  the  usual  physical  situation 
created  by  placing  a  len3  or  an  aperture  in  an  infinite  screen  in  the  field  of 
a  source.   In  such  a  problem  interest  centers  on  finding  the  field  transmitted 
through  the  lens,  that  is,  the  field  in  the  image  space. 

An  exact  solution  of  the  specific  lens  problem  or  of  the 
generalization  herein  considered  is  not  known.   The  problem  is  therefore  modi- 
fied somewhat  in  this  article  so  that  at  least  the  reformulated  problem  can  be 
attacked.  The  general  theory  of  asymptotic  evaluation  of  steady  state  fields 
is  then  applied  to  the  new  problem. 

3.2.V/e  define  an  optical  instrument  as  a  medium  which  is  homogeneous 
in  two  half  spaces  z  <  z  and  z  >  z.  and  which  may  have  arbitrarily  varying  but 
sectionally  continuous  £,  LX   ,  and  CT  in  the  domain  z  <z  <  z.  .  The  domain 

z  <z  <  z.  is  referred  to  as  the  lens.  Let  us  assume  that  £  ,  M  ,  C  have  the 

2  ? 

constant  values  €=  n     ,  f*  =1,<7~=  0  in  th3   "object  space"  and£=n  tf4=lt(7~>=  0 

in  the   "image  space". 

— >         -j» 
The  field  E  and  H  established  in 

the  image  space  by  a  periodic  source  at  P     will 
itself  be  periodic   (see  article  2)  and  will  satis- 
fy Maxwell's  equations. 


(3.1) 


curl  H —  E. —  E  =  0 

c   t     c 


curl  IS  +  ■*=-  E  =  0 


c    t  .  Figure  S 

in  source-free  regions.   If  we  write  E  and  H  as 

6.  See  the  reoort  mentioned  in  footnote  1. 
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Real  part   of  uU.y.sJs1      "  and  Real  part  of  f(x,y,z)e"   *°bl    respectively, 
where  u  and  v  are  complex  vectors,   then  u*  and  v*  are  themselves  solutions  of 

the  time-free  complex  Maxwell  equations 

2  --- 
curl  v  +  ikn  cj   =  0 

(3*?)  curltf-ik?-     =0 

We  assume  that  there  are  no  electric   sources   in  the  image  space.      It  is  assumed 

also  that  u*  and  v  are  regular  everywhere  in  the  image  space   z  >  z   . 

The  requirement  that  u  and  v5* are  regular  solutions  of  (3.2)  is  of 
considerable  physical  significance.  According  to  geometrical  optics  an  optical 
instrument  may  have  a  perfect  focal  point  at  a  point  P  of  the  image  space,  i.e. 
a  point  towards  which  all  rays  from  P  converge.  The  assumption  of  regularity 
precludes  a  singularity  even  at  ?  and  hence  implies  the  existence  of  a  limit  of 
resolution.  ITo  infinite  energy  densities  will  be  possible  at  any  point  of  the 
image  space. 

Let  the  time  periodic  source  operating  at  P     be  characterized  by 
F  =  g(x,y,z)e~  (see  article  2),   where  '£=■  0  with  exception  of  a  small  neigh- 


- 


borhood  of  P  .  We  consider  the  associated  "pulse  source"  F  =  g(x,y,z)^(t)  where 

£n(t)=0  for  t  <  0  and  lj(t)=l  for  t  >  0  and  assume  that  E  (x,y,z,t)  and  a"  (x,y,z,t) 

represent  the  electromagnetic  pulse  propagated  by  this  source.   This  field  will  be 

discontinuous  on  the  branches  of  a  hypersurface"^'(x,y,z,t)  =  0,  the  branches 

being  represented  by  t  =  l/c  "/j!(x,y,z).  In  the  image  space  where  £  ,  M,   and  (7~ 

are  continuous  (in  fact,  constant)  and  there  are  no  sources  it  follows  from  our 

first  paper  that  these  brandies  are  wave-fronts  and  so  must  satisfy  the  eiconal 

p  p  p  p 

equation*"^  +\4r    ■*■  "\£"  =  n  .  These  wave  fronts  characterize  the  internal  re- 
*  x  •   y   *  z 

flections  in  the  lens  of  the  optical  instrument. 

If  the  functions  €,  li,   ©"of  the  optical  instrument  are  sectionally 

analytic  and  if  the  source  is  given  by  a  sectionally  analytic  vector  field  F  then 

we  can  assume  that  the  pulse  solution  (see  article  2)  E  and  E  has  sectionally 

oo 

analytic  components.     More  specifically  let  us   assume  that  the  part  of  the  x,y,z,t 
space  which  belongs  to  the  image  space  z  >  z1    is  divided,   by  the  hyper  surfaces 
r&=  ct,    into  domains  in  whose   interior  E     and  H     possess  derivatives  of  any  order. 

Then  by   (2.3)  we  na^'  represent   the  solution  u  and  v  of  equations  (3.2)  by  an  in- 

V"  •*•    ik '  — '  -      i '-  *  :  '■ 

finite  series  of  Particular  solutions  of  the  form    JL.  IL-e  andx_V,~  e      '  ~~ 

respectively,  where  the  summation  over<s£  is  a  summation  over  the  many  wave-fronts 

-> 
passing  a  given  point  x,y,z  of  space.     As  remarked  in  article  2,    each  of  the  IL^ 

-^» 
and  Vcs  is  an  infinite  series  in  the  basic  variable  A  and  the  first   terms  in  the 

infinite  series  for  u"  and  v  are^p  |      e         "^  and<SL|li  J     e       ^':  ,    respectively. 

These  sums  approximate  u  and  v  for  small  wave-lengtns. 
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We  consider  in  particular  the  part  of  each  first  term,  namely 

"*•    ik"\£*         "■**   ik"sA 

Uqq  and  VQe  which  "belongs  to  the  directly  transmitted  wave- front 

*)^*(x,y,z)  =  ct,    since  it  alone  contributes  to  the   image  formation  of  the 
optical  instrument.     We  assume  that   these  two  quantities  are  known,   that   is 
we  assume  that  the  geometrical  optics   solution  of  our  diffraction  problem  is 


known.      This  means   that  we  know  ^3_ 


---- 


ik>^"   ,  rfri   ik-»A   ,   r-fi 

a^c.  _ii  wo      ,  where£2  J  and 


JE  j  are  the  electromagnetic  vectors  of  the  pulse  solution  on  the  extreme 

"boundary  of  penetration  of  the  transmitted  wave-front  >^(x,y,z)  =  ct.   Oar 
problem  thus  "becomes  to  find  a  solution  of  Maxwell's  equations  (3.1)  which 
in  case  k  ->  co  has  this  given  asymptotic  "behavior. 

The  "u  and  V  we  Seek  are  thus  far  not  sufficiently  limited  to 
determine  a  unique  solution.  On  the  other  hand  the  conditions  presented  "by 
the  presence  of  the  medium  between  z  and  z..  are  too  difficult  to  treat  exactly. 

We  therefore  rephrase  our  problem  and  thus  the  conditions  on  u*and  v*to  define 
a  more  amenable  problem.  The  significance  of  the  new  conditions  will  be  more 
apparent  later  on. 

3. 3. The  wave  f  ronts  y^=   ct  of  the  directly  transmitted  wave  are 
determined  by  the  rays  from  the  source  at  ?  .  These  rays  form  a  two  para- 
meter system  of  straight  lines  in  the  image  space,  and  the  wave  fronts 
~)r   =   const  are  normal  to  these  rays. 

We  now  consider  not  only  those  straight  line  sections  which  lie 
in  the  half  space  z  >  z.  but  their  total  extension.  This  means  that  we  supple- 
ment the  image  space  by  a  "virtual"  extension  to  a  total  homogeneous  space. 
The  rays  of  the  optical  instrument  cover  a  certain  part  of  the  total  image 
space  which  is  called  the  field  of  the  instrument.   In  particular  we  call  the 
field  regular  in  those  parts  where  there  is  exactly 
one  ray  through  each  point.  Obviously  the  field 
is  not  regular  in  the  neighborhood  of  focal 
points  or  caustics.   In  fact  two  or  more  rays 
go  through  the  same  point. 


Figure  9 
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We  shall  confine  our  discussion  to  systems  of  rays  which 
satisfy  the  following  conditions: 

1.  Che  ray  system  is  nowhere  telescopic  i.e.  not  more  than 
one  ray  shall  have  the  direction  p,q,r.  There  shall  he  no  "returning"  rays 
i.e.  r  2  0.  A  ray  thus  can  he  characterized  uniquely  hy  the  coordinates  p,q. 

2.  There  exists  a  finite  domain  of  the  x,y,z  space  outside  of 
which  the  rays  form  a  regular  field  and  the  associated  wave  fronts  satisfy 
the  conditions:  The  Gaussian  curvature,  K,  on  any  particular  wave  front  is 
positive:  X  >  0.  The  wave  fronts  are  convex  towards  the  origin  x=y=z=0  and 


KR  v>  1  if  E  =  a  <-y  vz    approaches  infinity  along  any  one  of  the  rays. 
This  means  the  wave  fronts  assume,  in  the  distant  parts  of  the  space,  more 
and  more  the  shape  of  spheres  concentric  to  the  origin. 

3.  There  exist  two  wave  fronts  "&~-   C.  and  "/r=  0     on  opposite 

sides  of  the  origin  which  can  he  embodied  into  a  closed  convex  surface, JL   , 
which  has  continuous  surface  normals 
and  whose  Gaussian  curvature 
K  I  ef   >  0.  (Fig. 10).  The  ex- 
terior normals  of  X3 determine,  outside  of 
— ."  a  regular  field  of  straight  lines 
which  contains  the  rays  of  the  instrument 
as  a  part.   Since  the  rays  are  perpendicular 
to  the  wave-fronts,  it  is  clear  that  any 
exterior  parallel  surface  JL,  of  JL  is  a 

closed  convex  surface  of  the  same  type.  Thus  an  infinite  set  of  closed  con- 
vex  surfaces  jn»  exists  which  contain  wave  fronts  of  the  instrument  as  parts. 
The  two  wave  fronts  which  comprise  part  of  X  will  he  denoted  hy  JS3'  and  12", 
and  TTJ  '  will  denote  that  section  which  lies  in  the  field  of  the  instrument. 

Consider  now  the  vector  fields  II  =  CUT  J  and  V  =  IH  ]  in  a 

o  o  o-o^ 

domain  of   the  x,y,  z  space  where  the  rays  of  the  instrument  form  a   regular 

field.     As  stated  ahove,    the  quantities 


Figure  10 


II  e       "^    ana  V^e 


0.3;  o  o 

determine  the  electromagnetic  vectors  ^and  v*  of  the  directly  transmitted  field 
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ia  the  case  of  small  wavo  length.  Prom  the  thoory  of  the  first  paper'  we  know 
that  in  our  homogeneous  image  space,  U  and  V  satisfy  the  conditions 


(3.*0 


=  o 

0       0 

grad)*v*  if  -  T  =  0 


and  the  differential  equations 

(3-5)  !   ,^. 

(_i  2  0 

along  the  rays.  The  parameter  s  measures  the  geometrical  length  of  the  rays.  By 
using  the  relation ^ refer  to  (5.22)  of  the  first  paper  and  use  the  facts "p*=  n  s* 
and'p  =  grad^i-^ 

(3.6)  ^^=-n     dl°fSl 

g 
where  K  is  the  Gaussian  curvature  of  the  wave  fronts  we  may  write  (3.5)  in  the 

form 

(3.7)  fc(,-=^-  fc/-==l-0 

and  conclude  that  the  vectors- 

(3.8)  -=r   ■  ^°»°->  J  ;=  =  &.*> 

are  constant  along  a  given  ray  p.q. 

3.^-.we  use  the  relations   (3.2)  for  the   definition  of  a  vector  field 

U  .   v      in  the  whole  image  space.      We  consider  the  exterior  normals   of  a  surface 

J?  as   defined  in  article  3.3.      In  the  exterior  of  the  surface  J^we  define  IT 

and  V     "by  (3.9)  if  a  point  x,y,z  lies  in  the  field  of  the  instrument  and  hy 

tj    =  v       =  0  if  not.      This  yields  a  sectionally  continuous  univalued  vector 
00 

field  in  the  exterior  of-^  . 

7.  See  also  Luneherg,   S.K. :  Mathematical  Theory  of  Optics,  Brown  University 
Lectures,   lghh,  Vol.    I,  pp.    25  and  H6. 

8.  In  the  first  paper  K  was  introduced  as  a  function  of  the  ray  parameter  C    , 
and  K(^»)  measured  the   expansion  of  a  wave-front  along  a  ray.      That  K(?"') 
can  he  identified  with  the  Gaussian  curvature  is   shown  "by  Silver,    S. : 
Microwave  Antenna  Theory  and  Design,  McGraw-Hill  Book  Co.,   1T.Y. ,  1°A9 »  P«lll+. 
formula  (22).      For  later  use  we  note  that   the  surfaces  involved  are  arbitrary 
provided  only  that   one  considers   expansion  along  a  curve   of  the  family  of 
curves  orthogonal  to   the  surfaces. 
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We  then  consider  the  interior  normals  of  J?  and  transfer  the 

vectors  of  the  exterior  field  along  these  normals  by  requiring  them  to  have 

the  values  given  "by  (3.7),  letting  A(p,q)  =  "L(p,q)  =  0  if  the  normal  does  not 

"belong  to  the  rays  of  the  instrument,  The  result  is  a  multivalued  vector  field 

.... 

u"  ,V  which  has  singularities  at  those  points  where  S  =  go,  that  is,  on  the 

caustic  of  the  ray  "bundle. 

Since  many  wave- fronts  may  pass  a  given  point  x,y,z  in  space,  in 

addition  to  the  transmitted  wave-front  thus  far  emphasised  we  obtain  more  than 

one  pair  of  vectors  u"   V  at  a  point  through  which  more  than  one  ray  passes. 

In  this  case  the  function  *^Kx,y,z)  is  multivalued.   If-^«,"^"1,  ...  are  the 

-">      -J*- 

values   of  "US  and  u '     V1    ;      u'",   V"     the  values  of  our  vector  field  ther 
»  o        o  00 


(3.9) 

0  0 

give  the  asymptotic  expression  for  the  solution  u^V  in  the  case  of  small 
wave  lengths  in  such  a  case. 

We  now  consider  Maxwell's  equations  (3.1)  and  (3.2)  not  only 
in  the  real  image  space  z  >  z.  but  in  the  total  homogeneous  x,y,z  space.  We 

wish  to  find  a  solution  u,v  which  in  case  of  large  values  of  k  approaches  asymp- 


totically a  given  vector  field  of  geometrical  optics  namely  the  field  (3.3) 
where  If  and  V  are  given  by  (3. S)  and  A  and  3  are  known  if  the  rays  form  a 


regular  field,  and  the  field  (3. 9)  if  not.  We  shall  require  in  addition  that 
u  and  v  are  regular  and  uniformly  bounded  in  the  whole  x,y,  z  space  and  con- 
struct a  solution  of  this  problem  in  what  is  to  follows. 

It  is  physically  plausible  that  the  approximate  solution  of 

— '  -  .-  — 
geometrical  optics  agrees  the  better  with  the  actual  steady  state  solution  u,v 

the  greater  the  distance  of  a  point  x,y,z  from  the  irregular  part  of  the  field 

of  the  rays.  Thus  we  are  led  to  the  conjecture  that  u*,v  attains  the  asymptotic 

— .^    ilr'V""'-*    ik"*^ 

form  U  e   ,  V  e     not  only  for  large  values  of  k  at  a  fixed  point  x,y,z 

but  also  for  large  values  of  x,y,z  and  a  fixed  value  of  k.   This  behavior  of 
u  and  V^is  indicated  already  by  dimensional  considerations.   The  degree  of 
approximation  must  depend  on  a  dimensionless  quantity;  we  shall  find  — — 

to  be  such  a  quantity.   Indeed  we  know  that  the  asymptotic  development  breaks 
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down  at  points  on  a  caustic  where  K  =   oo.     Vice  versa,    the  approximation  for 
a  given  value  of  k  will  he  best  for  smaller  values   of  K,    i.e.,    in  the  distant 
parts-  of  the  image   space,   for  K  — >- o  as   our  wave-fronts  become  spherical  at   large 
distances  from  the   source. 

Vfe  shall  use  the  above  hypothesis   for  the  establishment   of  the 
solution  of  our  problem  in  the  form  of  complex  integrals.     Afterwards  it  will 
be  shown  that  the  hypothesis  was  justified  by  proving  that  the   integrals  have 
the  desired  asymptotic  behavior  at  any  fixed  point  x,y,z  and  for  large  values 
of  k,   as  well  as  for  fixed  k  and  for  x,y,z  at   large  distances  from  the  source. 

3.5. The  problem  then  is  to   construct  a  regular  solution "u^v^of  the 
equations   {]>.?.)  which,   at  the  infinity  of  the  x,y,z  space  or  for  large  k  has  a 
given  asymptotic  behavior: 

*   -t   eik^0 
0 

->      £>    iky      „ 
v     -  v     e         ->  0 

0  ^ 

where  ^(x,y,z)  =  ct  is  the  equation  of  the  family  of  directly  transmitted 
wave-fronts.  By  using  the  relations  (3.8)  we  formulate  this  requirement  in 
a  more  concise  manner  as  follows:   If  00  is  approached  along  a  ray  p,q  of  the 
instrument  then  there  exists,  to  any  £  >  0,  a  sphere  of  radius  R     around  the 

Q 


origin  such  that 
(3.10) 


for  all  E.  >  fi 
0 


->    -ikl/-       . 

^ £(p,q)l  <    6 


If  infinity  is  approached  on  a  surface  normal  of  a  surface  jC* 
which  (normal)  does  not  lie  in  the  field  of  the  instrument  then 

R  |u  1  <  6. 

3,11  *\?l  <e 

if  H  >  R  . 
0 

finally  we  required  that  the  solution  u*,v  of  our  problem  satisfy 
a  "radiation  condition. " 

If  co  is  approached  along  a  ray  of  the  instrument  or  along  ary 
other  surface  normal  of  2?  then,  if  f  be  the  parameter  along  the  ray  or  normal, 
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(3-12) 

g 
or  explicitly^ 


(3.13) 


s  1    du"         ..    2    ->,!  <  c 

I  -&r  "    ika     u  *  -  c 

2  -:> !  >.     - 

-i    — -.—      -     Ikn    v  /  ^  c 


y  "  y 

These  last  conditions  are  generalizations  of  the  so-called  Sommerfeld  condi- 

tion  of  radiation.      They  insure  that  u,V represent  electromagnetic  waves 



which  have  the  same  propagations!  character  as  a  pulse  traveling  along  the 

wave-fronts"V"=  ct. 

We  shall   therefore  now  seek  a  solution  of  the  equations   (3.2) 
which  satisfies  the  conditions    (3. 10)   to  (3.I3). 

U.   Integral  Formulation  of  the   Diffracted  Field. 

h.l.We  now  proceed  to  solve  the  reformulated  diffraction  problem. 
3riefly,    the  accomplishments  of  this   article  are  as  follows.      In  the  homogeneous 
image  space  and  the  extension  of  it   hypothesised  in  the  preceding  article, 
Maxwell's  equations   reduce  to  the  wave   equation.      Green's  theorem  can  then  he 
employed  to  express  u  and  v"  as  integrals.     Further  analysis  transforms  the  inte- 
grands  so   that  they  involve  the  functions  A"(p,q)  and  B(p,q)  v/hich  are  supposed 
known  from  a  geometrical  optics   solution  to  the  problem.     The  final  integral 
formulation  is  given  by   (h.  15). 

Later  articles  will  evaluate  the   integrals  for  u  and  v  asymp- 
totically and  show  that  they  have  the   expected  behavior  for  large  k  or  for 
large  distances  from  the  source. 

9.      If  \J/  (x,y,z)  =  ct  is  equation  of  a  family  of  wave  fronts,   and  x(  ?),   y(  T), 

s(T)   the  equations   of  an  orthogonal   trajectory,    then  \j/    ,  \if     and  W     are 

2  2  2^2^  Z 

proportional  to  x',y!,   and  z1.      Since  "\j/     +  1/      +  1/       =  n       if  we  choose  f  so 

Q  O  O  O 

that  x     +  y     +  z     =  n',    then  W     =  x,   V/     =  y,    W    =  z.     Hence 

■  x  y  *  z 

du       _ 


— -    -  u    \y    +u    w     ■*•  u    \y 

dT  XII  y    V  y  z     r 


z 
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h.2.  To  derive  the  integral  representation  for  the  vector  fields 
u  and  v*  we  eliminate  u  from  the  Maxwell  equations  in  the  usual  way  (3.2)  by 
taking  the  curl  of  the  first  equation,  substituting  in  the  second,  expanding 
curl  curl  V*  as  -  A  v*"  +  W*   v*  and  using  the  fact  that  V  •  V=  0.   One  obtains 
the  wave  equation  for  v.  Likewise  one  obtains  the  wave  equation  for  u. 
Specifically, 


(Kl) 


2     2  ■*■ 
k     n    u     =0 


v  +  k2  n2^    -  0. 


We  now  apply  Green's   theorem  for  the  solution  of  equations 

(U.l).      Consider  the  equation  for  Tu(x,y,z)  and  a  domain  33  of  the  xj,z  space 

bounded  by  a  closed  surface  1.  .      We  remove  from  &  a  small   sphere  of  radius 

around  a  fixed  point  rjpZ,    of  D  and  apply  Green's   formula  in  the  remaining 

part  33  -6  to  the  function's  and  the  function 
inkR 


G  = 


E 


where  3.  ="/(x-x1)2+   (y-y^)2*   (z-z^2   . 


Green's   theorem  states  that 


(l+.2)y7/'(uA&-GAu)dxdydZ  =/ftfi  |i  -  G  |§)d<y. 

■13- £  nr 


S:  £ 


The  surface   integral  is  over  the  surfaces  X"  and 

Figure  11  € ;    the  derivatives     ™—     and  ——-    denote 

o'  -J  cJZs 

differentiation  in  the   direction  of  the  exterior  normal  of  the  boundary  of  D  -  £  ; 

and  d  co  is  an  element   of  surface. 

10 


It  follows  in  the  usual  manner   letting  £"^"0  that 


^.^^=#^|§  "G  |?)*^« 


(U.3) 

Since  u  is  regular  in  the  whole  space  we  conclude  that  (U.3)  holds  for  any  closed 
surface U  around  the  point  x.,y1,z.. 

We  can  choose  in  particular  any  exterior  parallel  surface  of 
^distance  -£  from  a  fixed  surface  J?n   of  the  type  assumed  in  article  3«  X*  as  well 

as  the  parallel  surface  X1  contains  two 
wave  fronts  V=  const.  Let  us  denote 
these  sections  by  X* ' »  X  "  andjC  ' .  3T  H 

respectively,  and  the  remaining  sections  by 
P'S'and  X'"  . 


10.  See,  for  example,  Slater,  J.C.  and  Frank,  K.H.:  Introduction  to  Theoretical 
Physics,  McGraw-Hill  Book  Co.,  F.  Y. ,  1933,  PP«  218-219. 
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related  hy  the  equation 


Since  corresponding  surface  elements   on^l*"    ancuL, 
11 


o     "  are 


KdCiJ  =  ::  - 

0     0 


We  can  replace  (H.3)  hy  the  integral  over  J^7  : 


0 


(U.U)  im?-  -  //  (£•      :  -  a  ;  -"-')   ^ 


In  this   integral  U,   G  and  E  are  the  values  assumed  at  the 

point  S  a£JL  which  lies  on  the  surface  normal  through  the  point  S 

o 

Since   (k.k)  is  true  for  all  parallel  surf aces  j2  it   remains 

true  in  the  limit  as^->co,    that  is,   we  have 


■  »      K 
(4.5)  Iffiu  =  -     lim      -'_■     (^-.-^_g      pi)     -2.  dca> 

"(0 


.  y  -  x   — 0 


Hence   the  limit   of  the  integral   (h. 5)  may  he  found  "by  determining  the  limit 
of  the  integrand  provided  that    the  convergence  is  uniform  on 


0 


inkH 
Now,      since  G  =     -2—-     ,  3    -  (lka  -     J)0    S    j 

and  since  the  wave  fronts  "become  more  soherical  in  shape  ,  ~—  '>1,     Also 

2  ®  * 

hy  an  assumption  of  article  3,   KB.  ~>  1     .  lurther,  the  integrand  in  (h.5) 

thus  has  the  same  limit  as  the  expression 

(h.6)  feikn-  h-  [   lr  e^' 


--.- 


Bat  "by  assumption  (3. 10),  since  K->  0   at  infinity,  u->0.  By  assumption  (3. 12), 
if  infinity  is  approached  along  a  ray  of  the  instrument 

(4.7)  a(§-  -  ikn^uJ-^O. 

12    a^>  _   2 

Since     — -^  =  n   we  may  write 


R  (dr  -  lk  or  *>  ->°- 


11.  See  footnote  S  of  article  3* 

12.  See  footnote  9. 
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tfe  may  change  parameters  from  T^to  1>   ,   the  derivative  in  the   direction  of  the 
normal,    and  cancel  the  factor  -~~     .     However,    (h.7)  mast  be  taken  in  the 

direction  of  a  particular  ray,      (see  the  assumption  ±i  article  3.3  that  the 

« — ,* 

direction  (p,q)  is  uniaue)  whereas  -    as  required  in  the  integral  (h.<j)  must 
he  outwardly  drawn  to*-  .  Eence  we  write 

CM)  s(      :  ik|||-k)*o. 

Hence   (4.6)  has  the  same  limit  as 


Ik    (ntlffllUZ*..*** 


(4.9) 

for  their  difference  approaches  zero. 

How     K<^|c    |gra<i^/j  =  a«      If  °^Q  chooses  a  ray  normal  to  the 
portion  ;; ' '   of  ^     along  which  to  approach  infinity  then  along  this  direction 

we  must   take  the  plus   siteft  in  (h.9),    for  we  have  denoted  "by""77'    that  section 

*  —  0 

of  y     which  lies  in  the  field  of  the  instrument  and  from  which  rays  of  this 
•*-  o 

field  emerge.      On  J  "    ,   however,   the  direction  of  approach  to  infinity  i3   the 
negative  of  the  exterior  normal's  direction  and  eo  (4.9)  vanishes   on  this  seg- 
ment of  J2      .     OnV"1   R  u  has  the  limit  0  "by  (3. 11)  and  so    (4,9)  vanishes 
there  also. 

It  follows  then  that  the  integral   (h.5)  need  "be  taken  over  I?' 

only  and  that  the  integrand  may  he  taken  to  he 

^•10)  2iknRu-K     eiknE         . 

0 

It  remains  to  determine  the  limit  of  En  e     .By  (3. 10)  and 

2 

the  assumption  that  KB  approaches  1,  we  have  that 

Bu--r(p,q)eik"^(SUo 
as  R  approaches  00  . 

Consider  nest  the  base  point  Q,  of  the  perpendicular  dropped  from 
(x-.y-.z.  )  onto  the  ray  through  S  (Figure  12).  We  have  the  relations  ~^{x,y,  z)=ct 
for  a  wave  front  and  n  =  —  where  v  Is  the  velocity  of  an  electromagnetic  wave 
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along  a  ray  in  the  medium  of  dielectric  constant  6.   Then  "*/^(x,yt  z)=nvt=n 
tines  the  distance  traversed  "by  a  point  following  an  orthogonal  trajectory 
to  the  family  of  wave  fronts.   In  other  word3"^(S)  is  the  optical  length  of 
the  path  from  the  source  to  S.  Hence 

V(S)  =  y^tQ)  *  n(SQ)  .X3 

3ut  SQ,  -  E-*-0,   hence  "^S)  -  ha  -5*YXQ.)l      that    is,    (U.10)  has  the  limit 
2i  knKX(p  Q)eik~^^^     U3  thus  0-btain  for  ^«5)  the  integral  formula 


and  similarly 

(H.u)  T--   &//*<,10.Ue*,l>zodc«4      . 

Since  Z  da»  =  Kd<-0  is  the  same  for  all  wave  fronts  it  is  clear  that  the 
0   0 

■  1 
integrals  (U.ll)  are  independent  of  the  particular  wave  front  i-  which  we 


have  chosen. 
© 

U.3«As  remarked  a'oove  the   function  "^(Q.)   determines  the  length  of  the 
optical  path  along  a  ray  p,q_measured  from  the  source  P     to   the  hase  point  Q,  of  the 

perpendicular  dropped  from 
x-i,y-i»z,    onto  the  ray.     Con- 
sider now  the  "base  point  Q, 

of  the  perpendicular  dropped 


P0  Source 

Figure  13 


from  the  origin  onto  the  ray. 


13.     The  quantity  n(SQ,)  is  the  optical  path  from  S  to  Q}   the  optical  path 
"being  in  general,   hy  definition,  J  n  ds,   the  integral  taken  along  a 


J  n  ds  =     a )  —  =  c  _/dt. 


ray.     Further 
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We  shall  use  the  fact   taken  from  the   theory  0f  optic3  that  the  length  of  the 
optical  path  from  Pq  to  Q,     is  given  "by  Hamilton's  nixed  characteristic  of  the 
optical  instrument.1^"     It  is  a  function  of  the  space  coordinates  of  the  point 
P     and  the  direction  numbers  ptq,r  of  the  ray  in  the  image  space  and  is  denoted 

W(?o,p,q)  =   V.'(p,q)  =V^(QQ)    . 

Now  since  *-,  -*»  and  —  are  the  direction  cosines  of  the  ray  p,q,r,  the  optical 

path  Q  Q  is  given  directly  "by  x^p  +  y,q  +  z.r. 

Hence 

(U.12)  ~Y"(Qj  =  VJ(p,q)  +  Xjp  +  yxq  +  Z]r  . 

With  the  aid  of  W  we  thus  may  express  the  u  and  v"  of  (U.ll)  in  the  form 

*<w)  =  _  to  ffti^c)  e^C^^yq^r) 
(U.13) 


2n 


2  d<^> 
0        o 


/     „     n  _       ikn    //  ^gv         n     ik(Y/+xp+yq+zr)  _ 
(x.y.z;  =  _  —  J^  B(p,q)  e  ^'H  Kq 


o 


wherein  we  have  dropped  the  subscripts  from  (x.  ,x  z.)   . 

How  the  surface  of  integration  can  "be  any  wave  front  of  the 
instrument  and  the  differential  Kd CO  is  independent  of  the  wave  front.     Hence 

we  may  determine  it  by  choosing  a  wave 
front  which  is  far  away.  Consider  two 
wave  fronts     p'and  T\     Me  have  the 

relation  Kdco=  K  d  co   and,   since,  by 

P.l.r  0  0* 

2. 

assumption,  KR-»-i 


K  d  co      =     lim 

0         °       -o 

R->-co 


•  ► 
d  <x> 


where 


:r 


Figure  lU 


R- 


VA  A     2 


Kow,   let  us  assume   that   the  wave  front  jT1     is  given  in  parametric    form 


(U.1U) 


£=    4    (p.q) 


lh.  See  Synge,  J.L.:  G-eometrical  Optics,  Cambridge  University  Press,  1937.  P.2U. 
or  Luneberg,  R.  K.:  Mathematical  Theory  of  Optics,  Brown  University  Lectures, 
p.  113. 
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that  is,  "by  a  vector  whose  components  are  functions  of  p  and  q.  She  wave 
front  X*  then  is  given  hy  a  vector 


n~    P 


S=       ^"(p8q)   + 
where  /is  the  constant  distance  of -X7  from  J1     and  p,   the  vector 

n  -  p  -  q       .      In  accordance  with  a  standard  result  of  differ- 


ential geometry  the  surface  element   d.co  of  jC  is  given  "by  the  differential 
do,  =   /£  x  £  I  dpdq     =     /  (  Z?"p  +     ^    p^)  X  (    zfq  *     ^    P^)/    dpdq   . 
Howeverp     1^~  ->1     If  B-2*oo  .     Hence,    it  follows   that 

A 

Z  dcu    =     lim         _-    =     _  /p    x  P    /  dpdq 
°       °       £->co       E^  n         p         q 


1_    dodo 

n 


V2  2       2 

*  n  -  p  -  q 

\!e  introduce  this   expression  in   (U.I3)  eud  ohtain  the  diffraction  integrals 
in  the  final  form 


->,            v              ik      ff            X(p,g)  ikO.Vpjri-qy+rz) 

uU.y.z)  =  -     -5=    JJ  '  '  '    e 


257     "D      }/~2       2       2 
vn  -  p  -  q 

(U.15) 


dpdq 


-g-/            x  _       ik    //              B(p9a)  ik(vr-*px+qy+rz )      .    , 

v(x,y,z)  =  -  —    /  /  r'     •    e  H,/  dpdq. 

^         D     -,/!> 2 2 
Vn  -  p  -  q 

The  domain  of  integration  D  is  the  range  of  p,q  values  which  "belong  to  points 

oa  T1    in  accordance  with  (H.lh).      How  J?i   is  the  wave-front  in  the  field  of 
o  0 

the  instrument  and  hence  is  determined  "by  the  aperture  of  the  instrument.      It 

would  consist  then  in  the  set  of  all  p,q  values  of  ray3  leaving  the  aperture 

of  the   instrument. 
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The  integrals  (1+.15)  are  the  basis  of  the  diffraction  theory 
of  optical  instruments.   In  the  special  case  of  a  perfectly  converging 

spnerical  wave  we  have  V  =  const.   In 
this  case  the  integrals  are  known  as 
Dsbya's  solution  of  the  wave  equation. 


U.h.  The  integral  representation  (h.15) 
for  ti.   and  v  was  arrived  at  by  using 
some  of  the  assumptions  made  in  the  pre- 
ceding article,  such  as,  for  exam-ole, 
Figure  15  the  regularity  of  u"  and  v"  in  the  whole 

space  and  the  asymptotic  form  at  infinity.   It  remains  to  verify  that  these 
vector  fields  actually  satisfy  these  conditions.   We  demonstrate  first  that 
it  and  t"  are  regular  solutions  of  Maxwell's  equations.  The  proof  that  IT  and 
v*"have  the  required  asymptotic  "behavior  at  infinity  and  "beyond  that  the 
originally  postulated  behavior  for  large  values  of  k  will  "be  given  in  the  next 
two  articles. 

These  integrals  are  no  an  exact  solution  of  the  diffraction  prob- 
lem.  The  assumptions  introduced  in  article  3,  which  simplified  the  problem 
somewhat,  have  led  to  this  result,  and  even  though  the  integrals  will  prove 
to  have  the  required  asymptotic  behavior,  one  must  not  conclude  that  the 
original  problem  was  fully  met.  However,  by  virtue  of  the  fact  that  the 
integrals  do  have  the  required  asymptotic  behavior  they  should  serve  at  least 
as  good  approximations  to  the  solution  of  the  original  -problem. 

It  should  be  noted  that  the  integrals  can  be  evaluated  directly  in 
only  a  few  simple  cases.  Hence  the  theory  which  follows  on  the  asymptotic  be- 
havior of  the  integrals  both  at  infinity  and  for  large  k  values  is  practically 
a  most  useful  approximate  method  of  evaluation.   As  remarked  elsewhere  the 
results  obtained  may  prove  useful  wherever  integrals  of  this  form  arise. 

The  proof  that  "d"  and  "v*  satisfy  Maxwell's  equations  (3.2)  is  made 
as  follows.  By  (3.U)  and  (3.8)  it  follows  that  A  and  xT  satisfy  the  condition 

p"x  3"+  n2r=  0 
(H.lb) 


p  x  r  -  f  =  0 


where  p0"  =  grad-^"  =  (p,q,r), 
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Let  us  now  form  curl  u"  and  carl  v"  by  differentiating  in  (h.  15) 
under  the  integral   sign,    remembering  that  A  and  3  are  independent   of  x,y,z. 
Then 

ik(V/*px+qy+rz ) 


curl  u  =  - 


(ik)2 


(U.17) 


curl  v  =  - 


2fT 


(ik)' 

2n 


P  a  A 


D     -j/2       2       2 
l/n  -  p  -  q 


// 


_2_. 


:' 


-,/~2      2      2" 

I'll  -  p  -  q 


dpdq 


^(W+px+qy+rz) 


Eence,    on  account  of  (h. 16), 


(U.1S) 


curl  v  +  ikn  u  =  0 


curl  u"  -  ik     v  =  0 


that  is,  Maxwell's   equations  are  satisfied  by  iT  and/T   It  is   clear  that  u"  and  T 
are  regular  in  the  whole  x,y,z  space.      Indeed  we  have  the  inequalities 


(U.19) 


1*1    < 


(v-/     < 


k_ 

2TT 


k_ 
2n 


Jf  /I(p,q)/ 

D      ,/~2       2        2 

Y  n  -  p  -  q 


dpdq 


qi*^^,, 


that  is,  u  and  v"are  uniformly  bounded  in  the  whole  x,y,z  space. 

5.  Asymptotic  Evaluation  of  the  Diffraction  Integrals  for  Large  k. 

5.1The  integrals  given  by  formulas  (h.15)  are  the  solution,  then, 
to  the  diffraction  problem  as  formulated  in  article  3.  Except  in  simple  cases 
it  is  not  possible  to  evaluate  them  directly.  Hence,  one  must  seek  other 
methods  and,  as  already  indicated,  this  paper  will  undertake  asymptotic  evalu- 
ation. 

To  obtain  the  asymptotic  expansion  of  these  diffraction  integrals 
for  large  k  values  it  is  necessary  to  relate  these  integrals  to  the  pulse  solu- 
tion discussed  in  article  2  and  thereby  to  the  general  theory  of  asymptotic  ex- 
pansions reviewed  there  and  more  fully  presented  in  the  first  paper.  This 
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relation  is  established  "by  means  of  a  key  step,  the  introduction  of  measure 

functions  ?  and  Q,  whose  derivatives  turn  out  to  "be  E  and  K  .   It  is  then 

—  °      ° 
necessary  to  obtain  the  discontinuities  in  S  and  H  . 

oo 

5. 2. The  vector  fields  u  and  v",   as  represented  "by  the   integrals 
(U.15),    can  he  interpreted  as  a  superposition  of  plane  waves   of  various 


directions,    amplitudes,    and  phases.      The  amplitudes  are  given  "by  the    ex- 
sessions  *  X(p^}  ^d  ^^     ' 


V~2        2        2 
vn  -  p  -  q 


V7T2 


;      the  phase  "by  "'(p.q),   and  the 


v   n  -  p  -  q 

directions  "by  the  portion  px  +  qy  +  rz  of  the  exponent.     For  convenience  we 
introduce  the  notation 

0U,y,z,   p,q)      5  px  +  qy  +  rz  +  W(p,q) 

For  the  asymptotic  development  of  u"  and  v"  the  integrals  (h.15) 
are  transformed  as  follows.  We  introduce  the  "measure  vectors" 


?(x,y,z,t)  =    ■— 


t< 


if . 

^ct        W  2       2        2 

v  v.  -  p   _  q 


(PjlSJ 


dpdq 


C5-1) 


£(x,y,z,t)  = 


2TTC 


Slp.q) 


OXcX,        i/~2~    2        2 
r  /n  -  p  -  q 


dpdq 


For  any  fixed  point  x,y,z  the  integrals  are  to  be  extended  over  the  part  of  the 
domain  D  of  p  and  q  values  where 

/        O  Q  /^ 

^  =  xp+yq+z    yn  -  p  -   q       +  W(p,q)    <   ct. 
If  x,y,z   is   a  point   of  the   real  image   space 
then, in  view   of   the  physical  meaning,^ 
is  greater  than  0. 

¥e   shall  now   shov;  that  the   inte- 


fitCt 


Figure  l6 


grals    (U.I5)   can  "be  written  as   sinrole 

15 

Stieltjes  integrals: 


u  =  -  ia> 


(5.2) 


-  iCOJ     e 

o 


I 


e     d?(t) 


i<»t 


dC(t). 


15.   It  should  perhaps  have  been  remarked  earlier  that  all  vector  equalities  are, 
a.s  usual,  no  more  than  3  separate  scalar  equalities,  one  for  each  the  com- 
ponents.  Thus  the  use  of  the  Stieltjes  integral  with  a  vector  integrand 
involves  nothing  new  as  compared  with  a  scalar  integrand. 


-■*!- 


The  identification  of    (5.2)  with.   {k.  15)  can  "be  made  as  follows, 
For  any  fixed  x,y,z,    the  function  <f>  depends  upon  p  and  q.      One  nay  then  con- 
sider the  region  D  of  p  and  q  values   to  "be  "broken  up   into  sub-regions  deter- 
mined "by  ct.    .      S     0<  ct.   where  t    ,    tn...t     ,,tn 
J        i-l  T  1  o'      1         n-1' 

are  increasing  values  of  t.and  t     and  t      (possibly 

*  0     n 

infinite)  are  the  extreme  values  of  t  which  need 

be  taken  for  a  fixed  x,y,z,  to  have  the  range  of 

p  and  q  values  over  which  0  is  defined  cover  the 

region  D. 

Figure  17 

The  choice  of  a  value  t  of  t  determines  by  (5.1)  the  function 

p(t4).  Let  us  consider  the  expression 

ik  0 


1 

(5.3) 


OV-ru^Je^i-i  , 


where  0  .    .    is  the  value  of  <p  determined  by  any  p.    . ,   q.    .    in  the   region  for 
which  ct     ..      j.   P  <  ct.    .     For  definiteness  p.    ,  ,    q.    .,    may  be   taken  so  that 
$  (P-    -it   °.-    -1  )  =  c"-    i    •     ^ou  *^e  expression  (5-3)i   ?part  from  a  constant 

factor,      is  the  value  of   one  term  in  the  sum  which  is   formed  to   obtain  IT  in 
(h.lR).      Indeed,    regarded  as  a  Stieltjes   integral,    the   integral    (U.  15)    is 


-»-  ^        ik     .  . 

u  =  -  2nc     —     lim 
2rr 

n-^-00       i=l 


[RV -*.„>].*•« 


=   -ico 


11.  tl^\)-^\^ 

n>oo     i-l 


ik  cti.x 


=  -icj/ 


Sj 


eiCOt     dPCt)     . 


The  integration  over  the  range  0  to  co    for  t  certainly  yields  the  domain  D  of 

(1+.15). 

The  identification  of  the  form  of  u"  and  v"  in  (5.?)  with  the 

diffraction  integrals  (U.15)  will  now  make  possible  the  application  of  the 
general  theory  of  asymptotic  expansions. 
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We  shall  assume  in  the  following  that  A(p,q),  E(p,q)  ana  W(p,q) 

are  integrahle  functions  which  are  sectionally  analytic.  Also  the  "boundary  of 

D  shall  consist  of  a  finite  num'oer  of  analytic  curve  sections.   It  is  clear  that 

then  P  and  Q,  will  he  continuous,  sectionally  analytic  functions  of  x,y,z,t.  We 

may  then  replace  dP(t)  hy  P.  dt.  Also  the  time  derivatives  P  and  0  are 

t  t  t; 

sectionally  analytic  hut  not   necessarily  continuous.      If  they  are   integrahle  we 
introduce 

]T(xf*,z,t)  =P^ 

(5.^)  Ho(:c,y,z,t)  =  Q^ 

and  may  write 


u  =  -   icu 


'       +      io»t    ,. 

E     e  at 

o 

0 


r   T-r       icot 
$  =  -  iw         He  dt 


(5-5) 


'0 

Formulas    (5.5)   should  now  he  compered  with  formulas    (2.1)  which 
relate  the  u  ana  v"  factors  of  the  steady  state  field  with  the  pulse  solution 

arising  from  the  same  source  distrihution.     Actually   (5.5)   i-s  &  special  case 

~*~  "*~ 

of   (2.1)  -orovided  E     and  H     are  0  at   t  =  oo    .      This   is   true  for  the  vectors 

o      o 

s"  snd  H  defined  hy  (5.4)  and  (5.1).   In  fact,  from  (5.1)  it  follows  that  for 

fixed  x,y,z,  P"  and  Q  are  constant  vectors  in  t  for  all  t  larger  than  a  fixed 
value  of  t  for  which  the  set  of  points  where  0  <  ct  fills  the  total  domain  D. 
Hence  3  =  P  =  0  anc  H~  =  Ql    =   0  for  sufficiently  large  t. 

-*.->-  ->.   -^_ 

It  is  also  true  that  P  and  Q,  and  hence  the  derivatives  E  =  P 

0    t 

and  h  =  cT  are  solutions  of  Maxwell  's  equations.   To  prove  this  we  introduce 
0     u 

the  functions 

-■2  2    -*- 

and  write   (5.1)   in  the   form  P  =     — -  P*   ;      <|"  =  ^—  Q*       where 

9  IT  dtd 


*"•£»%  ,    t(??q)  X(t-^)dndo, 

Vn  -  p  -  q 


*  ■    h  if         ?(p'q)        -V<»  -  i  «•  J  ***    ■ 


D  V  2        2        2 

Vn  -  p  -  q 


-33- 


That  P  =  P*     follows  by  differentiation  under  the   integral  sign  and  using  the 

definition  of     (t  -     -^  ). 

c 

We  show  first  that  P*  and  Q*  satisfy  Maxwell's   equations.      Since 

(t  -  — -  <p)  has   continuous   derivatives,    differentiations   of  P*  and  Ci*  v/ith 

respect  to  x,y,   and  s  can  "be  carried  cut  under  the  integral  sign.     Carrying  cut 
the  differentiations  called  for  in  curl  P*,   rearranging  terms,   and  using  the 
fact  that  grad  <p  =   (p,q,r)  =  p,  ue  obtain 


c  curl  P*  =  - 

2nc 


hi/     r*  r      **«  -  -«»  ** 


2        2       2 
a  -  p  -  q 


^"    "     ^   ^  /2       2        2  (  '     ^      ^ 

Vn - p - q 


->-  g  (5* 

and  two  similar  relations  for  c  curl  0*  and       ?* 

at 

Hence  by  (U.lo) 


(5.6) 


curl    "<$*     -     —    P*     =  0 
c         t 


curl     P*     +     -    Q*         _ 

C  Ti        =    0 


that   is,  Maxwell's   equations  hold  for  P*  and  Q*.      By  differentiating  these   equa- 
tions  twice  with  respect   to  t  we   obtain  the   same   equations   for  P  and  Q,  and  by  a 

third  differentiation  Maxwell's  equations   for  B     and  H     everyv/here  where   2*~   and 

oo  o 

-^ 

H     are  seetionally  analytic. 

We  thus  have  found  that  the    "measure  vectors"  P  and  Q,  and  their 

derivatives  E     =  P .  ,   H     =  Q,,    can  be   interpreted  as  pulse   solutions   of  Maxwell's 
o  t        o  t 

equations.      As  noted  in  article  2,    the  asymptotic  behavior  of  the   steady  state 
solution  "u*",v*" is    determined  by  the   discontinuities   of  the  associr-ted  pulse  solu- 
tion 3   ,   H   ,    that    is,    in  our  case  by   the  discontinuities   of  the    derivatives   of 
oo 

?",^1      We   obtain  the  desired  asymptotic  development    in  the  same  manner  as  before. 
i;e   introduce  the  vector  fields 

i     f00      dT(^) 


M(t)   = 
(5-7) 


^Jo  (7^-t)2 


2niy°        (r-  t)2 


_3U- 


As  pointed  out  in  article  2  these  vector  fields  are  analytic  functions  of  the 
complex  parameter  t  =  t.  +  it_  and  are  regular  in  the  whole  t  plane  with  ex- 
ception of  the  positive  real  axis  t.  =  0.  The  discontinuities  of  the  deriva- 
tives of  P  and  Q  determine  the  location  and  the  type  of  the  singular  points  of  V. 


and  h   on  the  positive  real  axis. 

By  application  of  the  theory  reviewed 
in  article  2,  a  representation  of  u**and  v"by  complex 
integrals  can  "be  obtained 


(5.s) 


u"=  J  Kit 

L 

v  =  J   K(t 


;  e     at 


J   e     dt  . 


U 


Figure  13 


The  points  t^  =  -"V^x.y.z)  where  P(t)  and  £(t)  have  discon- 
tinuous derivatives  deteraine  the  branch  points  of  the  Riemannian  surface  of 
the  vectors  M  and  h.     Again  in  accordance  v/ith  the  theory  of  article  2,  "by 
deformation  of  the  "branch  lines  of  this  surface  and  of  the  path  of  integra- 
tion  we  find  the  following  representation  of  u  and  v: 

u  =  Z_Uc<  e      ;  v  =  Z^V^  e   ro< 


(5.9) 


£< 


'•/here  U^    and  V^    are  given  "by  the   integrals 


dr 


(5.io) 


w 


.  /» 


(±y^.7-)eico'7ar 

-j    c 


-^ 


The  path  of  integration,  L,  is  now  a 
curve  of  the  complex  7~=   71  +  i  7^ 

plane  which  surrounds  the  positive 
imaginary  axis. 


■igure  13 
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We  shall  see  now  that  E  and  H  possess  discontina.it ies  of  the 

o  o 

type  which  will  permit  application  of  the  general   formula  (2.5). 

6.  Analysis   of  the  Discontinuities  of  ?!(t)  and  cf'(t). 

6.1.   The  method  we  shall  use  to  determine  the   singularities  of 
p'(t)   and  Q,'(t)   is  analogous   to  the  method  of  stationary  phase.      The  analysis 
depends  upon  stationary  values  of  the  function  <p  =  px  +  qy  +  rz  +  W(p,q)  where 

V/e  shall  therefore   investigate  the  function  0  somewhat. 


_V~2   2   2~ 

=  vn  -  v  -   q 


r  =  v  n '  -  p 

In  accordance  with  the  facts  introduced  in  article  h,    the 
geometric  meaning  of  $>  is  the  optical  length  of  the  ray  from  the  source  P 

to  the  base  point,  Q,,  of  the  perpendicular  dropped  from  x,y,z  onto  the  ray  p,q. 

Consider  now  a  fixed  point  x,y,z  which 

lies  in  the  regular  part  of  the  field  of 
(o.a.r 


''rays.  Let  vKx.y.z)  =  C  "be  the  wave 


front  through  x,y,z,  and  let  P  he  the 

intersection  of  the  ray  through  Q,  with 

~y^".  Since  "vK?)  J-s  tne  optical  path 

from  ?  to  P,  (p  =  C  ■*■  n  PQ.   It  is  geo- 
Pxgure  20  o     '  r 

metrically  evident  that  PQ,  and  hence  <fi  attains  a  stationary  value  for  a  ray 

which  passes  through  x,y,z.  Obviously  PQ,  =  0  and  thus  0=7^(x,y,z)  on  -phis 

particular  ray. 

We  may  thus  determine  the  ray  which  goes  through  x,y,z  by  the 

conditions  for  a  stationary  value,  namely, 


(6.1) 


0      =  x  -  2.  z  +  Y/   =0 
P      r      p 


0q  =  y  _  a  z  +  Wq 


0. 


If  x,y,z  lies  in  the  regular  field  of  the  rays  then  a  'unique 
solution  p(x,y,z),  q(x,y,z)  of  (6.1)  exists.  However,  the  equations  (6.1) 
are  true  also  in  the  irregular  part  of  the  field.   In  this  case  several 
directions  p(x,y,z),  q(x,y,z)  are  obtained,  that  is,  several  rays  pass  through 
the  point  x,y,z. 
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Consider  the  planes 
(6.2)  0(x,y,z;   p,q)     =  px  +  qy  +  rz  +  W(p,q)     =  C. 

These  planes  form  a  two  parametric   set  depending  upon  the  parameters  p,a. 

2  2  2  2 

Since  <p  ^  +  0     +  0z     =  n    we  recognize  that    0(x,y,z;  p,q)   is   a  complete 

integral  of  the  equation  0,  *  0>y+  0Z=  n2  and  hence  that    (6.2)   is  a  two  para- 
meter set   of  wave  fronts.     The  envelope  of  the  planes   (6.2)     thus  must  "be 
again  a  wave  front  y£-(x,y,z)  =  C.      We  obtain  yKx,y,z)  in  the  usual  manner 
for  determining  the   envelope  of  a  family  of  surfaces,   namely,   by  determining 
p  and  q  from  the  equations  (6.1)  and  "by  introducing  the  result  in  £)(x,y,z;  p,q). 

This   function y^(x,y,z),    of  course,    is  the  wave  function  of  cur 
system  of  rays.     The  set   of  planes    <p  (x,y,z;   p,q)  =  C     thus   represents  the  tan- 
gential planes  of  the  wave  front  y^(x,y,z)  =  C. 

We  determine  for  a  later  application  the  Gaussian  curvature  of 
the  wave  front  V""(x,y,z)  =  C.      We  know  that   this  surface  can  "be  found  "by 
eliminating  p  and  q  from  the  equations 

0p     =  °»       A     =  °»       &  =     C 
or  by   (6.1) 

x-2.z  +  W       =o 

r  p 

(6-3)  y-Is+v      =0 

r  q 

px  +  qy  +  rz  +  W  =  C. 
We  also  may  use  these  equations  to  determine  x,y,z  as  functions  of  p,q.  This 
gives  us  the  parameter  representation  of  the  wave  front ->^"=  C: 

x=-W  +  E-  (C  -  W  +  p  K  +  a  W  ) 

(6.U)  y  =  -  W  +  %  (C  -  W  +  p  W  +  q  W  ) 

q    nd  p     q 

z=         £-(C-W  +  ptf+qW). 

n2        q    q 

The  surface  element  dco   of  the  wave  front >^"=  C   is  given  by  the   expression 

d  cx>    =  /  X"  X   X"  I  dpdq 
P        q 
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where  X(p,q)  =  (x,y,z)  is  given  by  (6.3)  or  (6.10.  By  (6.3)  it  follows  that 


p 

X      -      "-    Z 

P         r     p 

+    <p 

PP 

=  0 

0 
y    -    *  z 

-.-     0 

=  0 

P    r  P      qq 


_  £  3  +  ^   =  0 

r  q     pq 


r  q    "qq 


the  quantity  0   ,  for  example,  "being  obtained  through  differentiation  of  (6.2) 

where  0  is  regarded  as  a  function  of  the  five  independent  variable  x,y,z,p,q, 

a     *        -/~2   2   2~ 
and  where  r=  '  n  -  p  -  q 


After  an  elementary  calculation  we  find  that 


(6.5) 


We  now  use  the  relation 


r  2 

dco    =     -  (    0         (2>       -    d)      )     dpdq   . 
n        rpp    V-qq         rpq'       -    H   ' 


pq 


n       r 


derived  in  article  \  and  conclude  that 

(6.6)  K  = 


r2(  0 


*pq 


pp  >"qq 
This  means:   If  the  wave  front  -y^x.y.z)  =  C  is  characterized  "by  the  equations 

0V     "  0,   0   =  0    ^(x,y,z;  p,q)  =  C 

then  its  Gaussian  curvature  K  can  he  obtained  with  the  aid  of  (6.6). 

6.2.  Consider  first  a  point  (x,y,z)  which  lies  in  the  regular  part 
of  the  field.  More  specifically  let  us  consider  points  at  the  left  side  of 
the  origin  so  that  the  wave  fronts  yr=   ct 
are  convex  towards  the  right.  We  saw 
that  at  a  fixed  point,  x,y,z,  in  the 
regular  part  of  the  field 

Cp=   px  +  qy  +  rz  •*•  W(p,q), 
regarded  as  a  function  of  p  and  q,  has 
exactly  one  stationary  point  in  D  and 
since  this  point  must  "be  a  minimum  point  it  follows  that  the  contour  lines 

$(p»o.)  =  constant  in  D  must  have  the  form  indicated  in  Figure  21.   These 

curves  surround  the  minimum  point  C   and  touch  the  boundary  of  D  at  at  least 

*  0 

two  -ooints  C,    and  C„.      We   find  the  point  C     "by  the   conditions    0      =    0       =0. 
12  0  P  q 

The  associated  minimum  value  of  ^  is   given  by  the  value  of  the  function  >^(x,y,z) 

at  the  point  x„y,z. 
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The  points  C.,  and  C  also  represent  stationary  values  of  0(p,q) 
provided  that  the  value  of  0  are  considered  on  the  "boundary  of  P.  Let  p  =  p(s), 
q  =  q(s)  "be  the  paranetric  representation  of  the  "boundary  of  P.   Then  on  the 
boundary  of  D  0(p,q)  =  £>(p(r),q(s)J  .  C^  and  C  may  "be  found  "by  the  conditi 


on 


(6.7) 


11  = 
ds 


rp  ds       q  ds 


The  associated  values  0  =  >^(x,y,s);  <fi=   >^:(x,y,z),  considered  as  functions 

?     2  ? 

of  x,y,z  must  "be  solutions  of  the  equation  -^-c  +  -y^  +  -^   =  n"  „-  indeed 

these  functions  are  defined  "by  C b. 7 )  as  envelopes  of  a  one  parameter  manifold 
of  special  solutions. 

The  vectors  P  and  Q,  defined  in  (5.1)  are  sectionally  analytic 
functions  of  t  if  we  assume  that  A(p,q)  is  analytic  and  that  the  boundcry  of  P 
is  en  analytic  curve.   Let  us  confine  ourselves  to  this  case.   The  only  singular 
points  of  P(t)  and  Q(t)  (i.e.  points  at  which  P  and  §"  are  not  analytic  though 
continuous)  will  "be  those  of  t  which  correspond  to  the  points  C  ,0.  ,C?  in  Figure 
21,  that  is,  the  values 

The  components  of  P  and  Q  thus  are 
functions  of  t  of  the  type  illustrated  in 
Figure  22.   Since  these  values  just  found 
are  the  only  ones  in  which  ?'  and  ~§}   are  dis- 
continuous our  general  result  (5.9)  allows  us 
to  conclude  that  u  and  v  must  have  the  form 

u  =  Ue  +  U.e  X+U?e         d- 


(6.3) 


=  Teik^+r   e11^^    eik^2 


where  U,  V,  IL.  ,  Vl  ,  u   v?  are  completely  determined  "by  the  values  of  P  and  Q,  in 

the  immediate  neighborhood  of  the  points  t  =  —  y^-,  t.  =  —  y^\,    t  =  —  ">^2  » 

respectively,  since  these  values  of  P  and  Q,  determine  the  singularities  of  M  and 
jf"and  therefore  of  the  U^  of  (5. 10). 


Ca 
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Figure  23 


CP'C  Consider  next  a  point  xsy,z   outside 

the  field  of  the  instrument.      Since  no   ray 

of  the  instrument  passes   through  the  point 

x8y,z  it  follows   that  the   equations  d>     =  <t>    =  0 

P    Q 

have  no  solution  inside  D,  that  is,  Q  (p,q)  has 
no  stationary  point  in  the  interior  of  D.   The 
contour  lines  0   =  const,  have  the  form  illus- 


trated in  Figure  23.   They  touch  the  "boundary  at  two  points  C.  and  C  in  which 

the  values  of  p    on  the  "boundary  assume 
extreme  values.   The  associated  values 

of  t  are  t1  =   J  V"i  and  X     =  i  V^  . 

The  components  of  P  and  Q,  have  singu- 
larities (i.e.  are  non  analytic)  only 
at  these  particular  values  of  t  (Figure  2h) 


and  hence  the  representation 


Figure  2h 


u"=  Ux  e      +  v*  e    c 


(6.9) 


v  =  V1  e 


ik>^l   2-  ik^2 


If  x,y,z  lies  to  the  right  of  the  origin  in  the  regular  part 
of  the  field  a  similar  pattern  of  the  contour  lines  as  shown  in  Figure  21  is 
ootained.   The  difference  is  that  C  is  a 

maximum  of  (p  .  The  components  of  P  and  Q, 
depend  on  t  as  shown  in  Figure  25;  u~  and  v" 
again  have  the  form  (6.8). 

Consider  finally  a  point  x,y,z 
in  the  irregular  part  of  the  field.   Since 
more  than  one  ray  passes  through  x,y,z  it 


follows  that  the  equations  t> 


$< 


Figure  25 
=  0  have  more  than  one  solution,  that 


p   ^q. 

is,  there  are  several  stationary  points  of  <0  in  the  interior  of  D.  These  stationary 
points  can  "be  minima,  maxima,  or  saddle  points  (Figure  26).  In  addition  to  interior 
stationary  points  we  obtain  a  number  of  stationary  "boundary  points.   The  totality 
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of  these  points  gives  the  set   of  points       t^    =     1  y£    where  f"(t)   and^t) 
are  not  analytic.      The  vectors  vTand  v"are  represented  "by  series 


a>    z  * 


J*Y* 


'ex. 


(6.10) 


of  finite  many  terns. 


_ 
n 


i 


ik  VS< 


In  what  is  to  follow  it  will  "be  assumed  that  x,y,z  lies  either 


outside  the  field  or  in  its  regular  part. 
Thus  we  do  not  discuss,  in  general,  the 
practically  important  question  of  the 
"behavior  of  u  and  v"  in  the  neighborhood 
of  a  focal  point. 

We  also  shall  not  abandon 
the  assumption  that  the  functions  A(p,q), 
3"(p,q)  are  analytic  and  that  the  "boundary 


I  c+ 


Figure  26 


of  D  is  an  analytic  curve.   It  is  clear  from  the  ahove  that  sharp  corners  of 
the  "boundary  produce  new  singularities  of  ?  and  Q~and  hence  additional  waves 


U 


ikVSr  .    ik^? 


<* 


*     "a 


in  the  representation  (6.S), 


6. 3. We  now  consider  the  singularities  of  p"(t)  and  §(t)  at  an  interior 


stationary  point.   Let  us  assume  that  <p  (p,q)  =  px  ■*•  qy  +  rz  +  W(p,q)  has  a 

O 


minimum  at  an  interior  point  p  ,q  of 

0      o 


^9o 


the  domain  D.     Let  ct     =    0  (p   ,q    )   = 

0      -^o^o 

-/■"(x.y.z)  "be  the  minimum  value  of  &  . 

The  vector  P(t)  is  given  "by  the  integral 

(5.1),  namely, 


Figure  27 


(6.11) 

the  equations 
(6.12) 


X(p,o) 


dpdq 


We  now  introduce,   instead  of  p,q,    the   coordinates   t,   9     "by 

P  =  /'U.9)  cos  0 
q  =    /^(t.e)  sin  0 


~;i- 


so  that   f  ,   0  are  polar  coordinates   around  the  point  p   ,q  .     A  contour  line 

<f>  =  ct  is  given  "by  its   equation  in  polar  coordinates/>-/©(t,Q).     For  t  -  t   >  0 

sufficiently  snail  a  closed  curve  is  obtained  which  is  described  once  if  9  goes 
from  0  to  2n.  The  components  of  the  vector  field  —  i£(p,q)  "become  functions  of 
t  and  0   , 

^Atp.oJ  =  r*(t,Q)  =  A*(^cos  0,    /osin  0) 

In  the  neighborhood  of   the  integral^ (6.11),   after  a  change  of  variable  which 

requires  the  introduction  of  the  Jacobian  of  the  transformat ion  .(6.12)  "becomes 

r  t   r2n 

(6.13)  P(t)  =    -~        I       r*(t,©)^      dtdo  . 

J  t      o 

0 

We  determine  the  function^(t,0)  from  the  equation 

ct  =  tfKp.q)  =   ft  ( />cos  0,^sin  0). 

We  may  expand 0(p,q)  in  a  Taylor  series  in/?-  f    since  p  and  q  are  functions 

of  f  .      For  example  ^  -   0p  P/0    +    0q  q^,  ^=    0?  p  pf>   +  (£pp  p2    +  ^  qpPf+ 

Also/°(t   ,0)  =  0  for^(t,<?)  is  the  polar  coordinate  equation  of  0=  ct  around  t   . 
o  '  o 

Finally,     <z5     =    <p     =Oatt=t     for  <t>    is  stationary  there.      We  have  therefore 
p  q  o 

in  the  neighborhood  of  t  =  t   : 

(6.1*0  c(t-to)  =  ^2  ^2(e)  *  />3    jz$(o)  +  ... 

where  0?(©),  ^,(e)  »•••  are  homogeneous  polynomials  of  order  2,3,..;'/)cos  0  and 
sin  0.   These  polynomials  satisfy  the  condition 

(6.15)  $4,(0 +  tt)  =  (-i)1'^2/(o)  . 

We  introduce  the  variable  s  =  c(t-t  )  and  consider  P   as  function 
/?(s,0)  of  s  and  ©  .  £y  (6.1*0  it  follows  that 

s  =  p  §  (/>,©) 

where 

^(/\©)  =  V02(O)  +^01©;  *  ... 

is  a  regular  analytical   function  of  f>  provided  that 

0o(O)  =4(0       cos2©  +  2    ^        cos  ©  sin  ©  +    &       sin2©)      >     0. 
r2  2     rpp  rpq  r  qq 
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This  however  follows  from  the  fact  that  &  -  (A  >  0  at  a  minimum 

'pp  ^q.q     ^pq 

point   of    ^(p,q.)  . 

By  (6.15)  we  conclude  that  j£~(  ^  s  Q  )  satisfies   the   identity 
(6.16)  #(-/>,©  *vr)     =<g{/o^)     . 

We  also   conclude  that  the  equation     s  =  /°jf(ptQ)  can  "be  in- 
verted in  a  unique  manner  so  that 

^=     oR(s,0) 
where  R  is  an  analytical  function  of  s. 

The  function  R( s,0)  satisfies  an  identity  similar  to   (6.16), 
namely,   R(-s,6  +  Tr)  =  R(s,0).      In  fact  from  s  =  />^(  />80)  it  follows  that 
(-s)  =   {- 1°)    <£  C~/>»e  *  Tr)  and,  "by  the  uniqueness   of  the  inversion, 
R(s,0)     =  R(-s,Q  +  TT). 

We  now  introduce  the  following  notation.     A  function  f(s,0) 
which  depends  only  on  even  powers  of  s  and  satisfies  the  condition  f(s,0  +tt)=   f(s,0) 
shall  he  denoted  hy  (G  ]    ;   a  function  f(s,0)  depending  again  only  on  even  powers 
of  s  but   satisfying  the  condition  f(s,9  +  TT)  =  -  f(s,0)  shall  be  denoted  by  /G,  J . 

Obviously:    [Gj  [gJ    =  [Gq7    ;    [&J  [dj  =    [c-J;  [gJ/gJ    =  [^  J  . 

With  this  notation  we   can  write 


and  conclude  readily 


f    =    S    [GrJ       *    UJ 

and  hence 

(6.17)  P/oQ  -  3  [Go]   +  [gJ. 

Consider  next  the   function 

~£*   (s,0)  =  jf*  (sB(s.O)  cos  08    sR(s,9)  sin  0)  . 

We  conclude  from  R(-s,0  «•  tt)  «=  R(s,G)   that     T^Cs.O  +  TT)  =     JC*(s,0) 

(6.1S)  T*  (s,0)     =  [G01  +  s  [Gj)   . 


By  confining  (6.1J)  and  (6. IS)  we  find 

!s  (8,o)/>  f3   =  S  [c-o]  +  [a1]  . 


/ 


\le  introduce  this  result  into  the  integral  (6.I3);      since 
2n  the 

[  S.j  dO    =  0  "by  the  property  of  G,    function,   we  obtain 


^t}    "    3Re/     *[(io)  ds    ~\  •*  +  <>? 


'0 


t  >*     , 


that  is,  P(t)  depends   only  on  even  powers  of  s  =  Yc(t-t    )     .      This  means  ?"(t) 
is  regular  in  the"  neighborhood  t  Z  tQ  of  t^o     A  similar  result,    of  course, will 
he  found  for  the  vector  ^ 

In  view  of  the  regularity  of  P(t)  for  tit     we  may  expand  it 
and  its   derivative  E     (also  Q  and  E   ;   in  a  Taylor's   series,   namely, 

K  •£  -h  LV»  (*  -  V 

(6.X9) 

where  [  B     J  and  [e   Jare  the  right  hand  derivatives  of  E     and  E     at  the  point 
o  ■*  0  00 

t  =  t   ,    that  is,   the  derivatives  ohtained  by  taking  limits  as   t  approaches  t 
o  0 

from  above.     Let  us  determine  explicitly  the  coefficients/"  B    7     and/  E~  ]       that 

__  °  0      °  0 

is,  the  discontinuities  of  the  vector  E  and  E  at  t  =  t  .  For  very  small 

000 

values   of  t  near  t    ,     A(p,q)  does  not  vary  much;      hence  we  may  say  fron  (6.11 ) 

0 
Since  the  Taylor's  expansion  of  <p  in  the  neighborhood  of  t  =   t     or   (p   ,q    ) 
begins  with  second  derivatives  of  0  with  respect   to  p  and  q  the  integral  on 
the  right  side  is  the  area  of  the  ellipse 

<PVV    (P-P0)2  *   2    0pq    (P-?0)(q-q0)   +     0qq(q-qQ)2  -    &><t-t0)    . 
It  follows  that 


ro^<*pp   4>m-& 


16.  P(t)  is,    of  course,    zero  for  t   <  t     . 


-Vi- 
and hence  with,  the  aid  of  the  result   (6.10)    : 

[Z0]    =A(p0,q0)  VT 
(6.20) 

iro]    =Blpo,qo)V^ 

where  K  is  the  Gaussian  curvature  of  the  wave  front  V"51  c*     through  the  point 
x,y,z   in  question. 

The  preceding  result  psrmist  us  to  derive  the  vector  fields  U,v 
in  (6.S)   .      We  apply  the  general   formula  (2.6)  to  our  case,    letting  m  -  X  =  0, 

D     =  [s  ]      and  C     =  0.      It  follows  that 
n       L  oJ  n  n 


y.  27  Sk  ;    v  -  I  [t"^ 


The  first  terms  in  these  series  are,  in  view  of  (6.20)  , 

(6.22)  ^Tjo    =Kpo,qo)YT 

7aBPJo   =^p0'%)^ 

We  must  now  recall  that  if  we  fix  on  any  point  x,y,z  in  a  regular 
part   of  the  extended  image  space,    there  is  one  ray  p    ,q     through  x,y,z  for  which 
(p  (p,q)   is  stationary.     Also  X(p   ,q    )  V~K  ,  which  is   constant  along  this  ray, 

gives   the  geometrical  optics  "behavior   of  instrument.     Hence  the  result:      If  PQ.qo 

is   an  interior  stationary  point   of  0=  px  +  yj  +  rz  •*•  W(p,q)   for  a  given  x,y,z 

then  the  associated     vector  fields  U   e ,   V  e        "    of  the   diffraction  integrals 

approach,   with  increasing  k,   asymptotically  the   electromagnetic   field  of  geometric 
optics,    namely, 

o  e  =     A^0'qo)l/K     e 

(6.23) 

=r     ikV_    -37  v    r=r      iky>- 

Vq   e  -     3(po,qo)VK     e 

Since  a  unique  p  ,  q  exists  and  passes  through  any  given  point  x,y,z  in  the  regu- 
lar part  of  the  field  (extended  image  space)  the  result  applies  for  all  rays  in 
the  regular  part  of  the  field. 

The  asymptotic  "behavior  of  the  two  other  vector  fields  in  (6.8)  shall 
"be  determined  in  the  next  section.  The  result  will  "be  that  they  approach  zero  if 
k  — »  co  . 
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6.U.V/e  now  consider  the  singularities   of  P~(t)  and  Q(t)   at  stationary 
points  of  the  "boundary  values  of  0(x,y,z;  p,q).      Consider  the   case  that 

<z*  =  px  +  qy  +  rz  +  W(p,q)  assumes  a  minimum  value  at  an  interior  x>oint  C     of  D 

o  * 

the  domain  of  integration  for  the  diffraction  integrals,   and  that   its  "boundary 

values  on   r  ,   the  "boundary  of  D,   have  a  minimum  at  C.    and  a  maximum  at  C 

(Figures   21  and  28).      The  contour  lines  <fi  =  ct,    and  <p  -  ct       touch  the  "boundary 

y  at   the  points  C1  and  Z^  .     Let  us   assume  that^2  +  rt2     f  0  at   these  points, 

so  that  none  of  the  points  C.    and  Cp   is  a  real   stationary  point  of  A  with  respect 
to  values  of  (p  in  the  interior  of  D.      Let  us  also   assume  that  the  curvature  at 
C^  and  C_  of  the  tangential  contour  lines  is     different  from  the    curvature  of 
the  "boundary  curve  at  these  points. 

The  vector  fields  £"(p,q)  andlT(p,q)   shall  "be  analytic  functions 
of  p  and  q  defined  in  a  domain  D1    somewhat  larger   than  D. 

In  the  neighborhood  of   t  =  t. 


'1  ' 
dpdq  - 


(6.23a)        fit)  =  gL  //  I 

The  first  integral  is  extended  over 
the  domain  in  D'  "bounded  "by  the 
closed  contour  line  fi  =   ct,  the 
second  integral  over  the  sickle 
shaped  domain  D, .   The  first 
integral  is  an  analytic  function 
of  t,  the  second  integral  not 
"because  of  the  sharp  corners  in  the 
boundary. *  The  singularities  of  P(t) 
thus  are  determined  "by  the  second  integral. 
In  the  neighborhood  of  t  =  t, 
respect  to  values  on  f  , 


1 
2nc 


//  I 


D, 


<ipdq  . 


Figure  28 
since  ft   has  a  maximum  there  with 


(6.23b) 


**>  ■  &  f{%  **  -  ik  f{  F 


dxido 


The  first   integral   is   extended  over  the    total   domain  D  and  thus   is   a  constant. 
It  follows  that  the   second  integral,    which  extends    over  the  sickle   shaped  Dp, 


Since,  for  t  <  t.  ,  the  integral  over  D,  is  0,  and  for  t  >  t,  ,  the  integral 
over  D.  is  a  non-constant  function  of  t,  the  left-hand  derivative  of  some 
order  of  this  integral  must  be  different  from  the  right-hand  derivative  at 
t  =  tr 
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*»-  fi 


Figure  29,   gives  the  singularities   of 
P(t)   in  the  neighborhood  of  t  =  %     . 

The  type   of  singularity  obvi- 
ously will  he  the   sane  at  C,    and  C      . 
Thus  we  may  confine  our  discussion  to 
the  neighborhood  of  t  =  t?  . 

We  choose  the   coordinate 
system  such  that  p=q=QatC_  and 

that   the  boundary  curve  T  is  tangential 
to  the  p-axis    (Figure  29).      Since 


Figure  29 


i  = 


QQ 


^  =  0   at  p  =  0,    q  =  0 


we  may  write 
(6.24) 


-        °      2    . 

q  -    —  p    r  q-jP- 


as  the   equation  of  the  boundary  curve:     K     is  its   curvature  at  C 

o  2 

The  equation  q  =  q(p,t)  of  any  contour  line  0  =  ct  may  be  obtained 

from  solving  ^=  ct,  with  t  constant  on  the  contour  line,  f or  q  in  terms  of  p.  We 

may  then  consider  the  identity  <23  (p,q(p,t))  =  ct.   For  fixed  t.   %^--d)  +  0       —L  =   0 

•   dp    ^p  'q   dp 


since  t  is  constant.  Then  since 


dp 


0  at  C     by  virtue  of  the  tangency  of  y  and  the 


contour  line,   it  follows   that     <p  (   o,o)  =  0  .      Since  we  have  assumed  that 

p 

2      2 
0    +  0       f     0  at  any  stationary  ooint  of  D,  we  have  that  0   (o.o)  4   0  . 
p      q  *q 

By  expanding  ct  =  0(p,q)  in  a  power  series  in  p  and  q  about  t  =  t, 
or  p  =  0,  q  =  0,  we  obtain  in  the  neighborhood  of  this  point: 


(6.25) 


c(t-t2) 


+  ttt  (  iO2  *  2^pq  ■*■  0nn<?)  *  .... 


2-H 


pp 


pq 


qq- 


q  rq 

The  contour  line  which  touches  the  point  C~  is  found  from  (6.25)  by  letting  t  =   t_; 


it  has  the  form 

(6.26) 

so  that   K  =  - 


0VV     2   . 


v> 


*-*-     is   its  curvature;      our  above  assumption  on  the   curvature   of 

0^ 


the  contour  lines  means  K  7*  Zn,   i.e.,  form  (6.24)  and  (6.26)  K^  +  ^JPP,    f     0. 


% 
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How  consider  any  contour  line  passing  near  "but  not  necessarily 
through  C_  .  V,Te  determine  the  abscissae  p.  and  p  of  the  points  of  intersection 
of  this  contour  line  ft  =  ct  with  the  "boundary  f* .  By  inserting  (6. 24)  into  (6.25) 
we  get 

(6.27)  (/!  30  (t~V  =  p2(l  *  \v   *  a2p2  +  "^  ' 

Since  the   right  side  is  positive  for  p  small,   an  intersection  is  possible  only 
for  such  values  of  t  -  t     which  make  the  left  side  positive.      Since  C0  is  a  maxi- 
mum of  0  on  the  boundary  f  we  have      0     <  0  near  t2  and  hence  t   <  t?  if  K     >K 
and  t    >  t     if  X     <  X.      This  result  is  geometrically  obvious.      We  next  introduce 
the  quantity 

(6.28)  s2    -    lfj|LE7    (t  -  ,  ) 

wq      0 
and  use  it   in  connection  with  formula  (6.27)  to  calculate  the  abscissa  v     and  v 

from  the  equation 

1  s  =  p  V(l   +  ^p  +  &2p2  f  ...) 

It  follows,  since  theoretically  a  solution  for  p  in  terms  of  s  exists,  that 

P2  =  s  R(s) 
(6.2S)  p1  =  -s  R(-b) 

v/here  R(s)  =  R  +  R^  +  ...  is  an  analytic  function  of  s.  p.  and  pQ  then  are  the  p 
values  of  the  point  of  intersection  of  any  contour  line  with  y». 


X     are  given  by  the  integral 


As  noted  in  connection  with  (6.23b),  the  singularities  of  p(t)  at 


(6. 30)      p>>  •-  ■&  JJ  £«*--  &{'»L  ^-J*-' 

2  "/0  X 

Now,    the   equations   (6.25)  and  (6.2S)   show  that   q(s,p)     is  a  function 
of  s   ;    it   follows  that   A(p,s),   which  arises  from  A(p,q(s,p)),    and 

r(p,s)  =    Vn  -  p  -  q  (s,p)     are  also  functions  of   s     .     ~^-  thus  is     s  times  a 

2  ° 

function  of  s   .     The  second  integral  in  (6.3O)   then  has  the  form 

s     rP2 

P2(t)  =  ■*■  ~    J       SdS       (  ^(p)  +  <3j(p)s2  fr^iv)^   +  ...)  dp. 

o      Jp_ 


17.  The  second  integral  may  be  obtained  from  the  first  by  the  change  of  variables 
P  =  Pt  1  =  q(3)«   ^e  Jacobian  of  the  transformation  is  then  '3Q-     . 


-us- 

We  conclude  from  (6.29)   that  the  integral  j       ^(p)dp     must  "be  an  odd  function  of 

Pl 

s  and  hence,    that  F~At)  has   the  form 

s 
P2(t)  =  +  J       s2   (<T   +  c]_s2  +  ...)  ds  =  s3(C^  +  C|s2  +  c|sU  +  ...) 


3y  reintroducing  -.he  variable  t  by  means  of  (6.28)  it   is  found  that 

3  ^ 

>2(t)  =  +vtpv    (r  +  Cj^ctg-t)  +  c2(t2 


(6.31)  I"  (t)  =  +vt^t    (r  +  3i(t_-t)  +  c"_(t„-t)2  * ...)  . 


This  gives  for  the  singular  parts   of  B     =  Pt(t)  and  similarly  of 

H     =  Q,    expressions   of  the  tyoe: 
0         z  7^" 

EQ  =  +  Vt2-t    (2_  —  U2-0"  )  +  Reg.   functions 

(6.32)  ^' 

]T   =  +  Vt2-t    (Ej£-  (t2-t)^   )  +  Reg.    functions. 

Let  us  determine  the  first  coefficients  C     and  D*    of  these  ex- 

0      0 

pressions.  From  (6.25)  and  (6.28)  we  get 

K"Ko  2 
q  =  -  — s-  s   ■»■...  .  From  (6.29)  we  have  p=  -  s  +  and 

P2  =  8  +  ...   .   Hence,  since  Pg(t)  =  ■—    f    '   ds  f   2  p -^  dp,  and 

Jo  J   P;l  r  a  s 

r 

since  in  the  neighborhood  of  U_  we  may  regard  —  as  constant  and  hence  remove  it  from 

c.  r 

the  integral, 

_  K-K         f   rs  fz  K-K       -r    , 

*>>  =  -    ST     R    ■«•/      *♦...--    js*  f .'    ♦.... 

-S  •' 

From  (6.28)  we  then  obtain 

^  K-K        r/2c(t-tp)       \^2 


ZK  3rrc     r    V   ^    (K  -K) 

^  q     0 

In  a  coordinate  system  not   specialized  as  in  Figure   29  one  has 


0      **"    0  .     Hence  it   follows  by  differentiation  of  pl(t) 

(and  similarly  of  Q"_(t))and  by  reference  to    (6.32)   that 


a-         1  a     r   2c  1 

tt  r    V  [K-K  I        ,  3/2        • 

0  (^P2  +  *\ ) 
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2     2 
Oar  result  is  not  valid  if  0     +  <fi       =0  i.e.  if  C  is  a  real 

stationary  point  of  <p  located  on  the  boundary  f  ;  it  also  breaks  down  if  K  =  K 

o 
i.e.  if  the  tangential  contour  line  p  =   ctp  and  the  boundary  curve  f  have,  at 

C~,  the  same  curvature.   These  exceptional  cases  require  a  special  consideration. 

We  use  the  above  results  to  derive  the  remaining  vector  fields 

^1'  %.'   ^2'  V2  in  ^.8).  We  know  from  (6.32)  that  at  t  =  t? 

(6.310  ^  =yVt~X"c^/  (t2-t)"   for  t  <  t2 

=  0  for  t  >  t   . 

The  zero  value  of  E  for  t  >  \,       holds  "because  t  is  a  maximum  relative  to 

values  of  0   on  f  and  so  is  all  the  more  a  maximum  for  0  for  values  assumed  by 

all  p,q  in  2.  (See  Figure  29  ).  Hence  P(t)  =    ff        -  dpdq  is  constant  for 

^     -^  <£  <ct   n  ' 

t  >  t_  and  E  (t)  =  P   =  0  for  t  >  t_  . 

<-  0         o  c 

In  the  neighborhood  of  t  =  t^  P(t)  is  given  by  (6.23a).  The 
first  of  these  integrals  is  analytic  and,  since  C,  is  a  minimum  point  for/>=  ct 
with  respect  to  values  on  T  ,  exists  for  t  <  t.  and  t  >  t..   The  integral  will 
be  denoted  by  1T( t-t^.   The  second  integral  is  sero  for  t  <  t.  and,  by  considera- 
tions similar  to  what  were  employed  to  deduce  (6.31)  and  (6. 3?)  gives  rise  to  an 
expression  for  EQ(t)  of  the  form  of  (b.32).  Hence  in  the  neighborhood  of  t  =  t.. 

it  =  SU-^)  for  t  <  t 

=  EU-tj)  -  Vt=t^L%^  (t-txf    for  t  >  tx 

Also,  the  function  R(t-t..  )  is  regular  at  t  =  t„  . 

We  now  apply  the  general  formula  (2.6).  The  result  is 
**"         C 

(6.35) 
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Ihe  main  terms  of  these   series  are 


"5 


(1)  -*<2) 


u 


1  2  (-ico)V2  2  2  (icofTZ 

where  C"^   and  ~CA   are  given  "by  (6.33)'  Tile  formulae  (6.35)  demonstrate  that 

U^,  Up  and  similarly  fi,  T    approach  zero  if  co  ->  co  ;  i.e.,  if  X  ->  0. 

We  know  <f>  (x,y,z;p,q)  has  exactly  one  stationary  value  in  D  if  x,y,z  lies 
in  the  regular  part  of  the  field  of  rays.  This  value  is  either  a  maximum  or  a 
minimum.   We  have  found  that  the  part  of  the  diffraction  integrals  (k.lo)   which 
is  characterized  hy  the  stationary  point  approaches  asymptotically  the  given  small 
wave  length  field  AVX  e  '    ,  B"yZ  e   'if  k  ->  co  .  The  remaining  parts  of 
the  diffraction  integral  converge  to  zero.  It  follows  that  the  solut ions u"",v~ given 
"by  (U.  15)  have  the  asymptotical  "behavior  which  we  have  postulated  originally. 

In  the  irregular  part  of  the  field  several  interior  stationary  points 
of  0  and  more  than  2  stationary  points  of  the  "boundary  values  may  "be  obtained. 
The  contributions  from  these  latter  points  converge  to  zero  if  k  ->  00  .   On  the 
other  hand  the  contributions  from  the  interior  stationary  points  give,  for  large 
values  of  k,  the  asymptotic  fields 

r*r«vsr   eik^'    *  jt*/^    eik^"'  +  ... 

r«fitf?     eikV"    +    fny-F    eik^"+... 

i.e.,  the  asymptotic  "behavior  required  in  (3.9).  This  is  a  consequence  of  the 
result  of  article  0.3  in  the  case  of  minimum  points.   The  same  consideration,  of 
course,  can  "be  made  in  case  of  maximum  points.   Saddle  points  require  a  slight 
modification  of  the  method  (Pt(t)  and  0^(0  have  logarithmic  singularities)  "but 
the  asymptotic  result  is  the  same  as  before. 

We  also  mention  the  fact  that  our  assumption  aoout  the  analyticity 
of  A(p,q)  is  not  essential.   If  jT(p,q)  is  sectionally  analytic  "but  discontinuous 

on  certain  curves  in  the  domain  2,  then 

these  curves  have  the  same  effect  as  the 

"boundary  curve  T  of  D,  that  is,  they 

introduce  new  points  whe  ."e  the  contour 

the 

lines  d>  =  ct  are  tangential  to/houndery 
curve. 

Figure  30 
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7.  Asymptotic  Behavior  of  the  Diffraction  Integrals  at  Infinity. 

7.1  Articles  5  and  6  analyzed  the  asymptotic  behavior  of  the  diffrac- 
tion integrals  (h.15)  for  large  k  or  small  A.  It  remains  to  "be  shown  that  the 
integrals  (h.15)  satisfy  the  "boundary  conditions  at  infinity  as  formulated  in 
article  J>.h.     We  do  this  simply  "by 
transforming  this  problem  in  such 
a  manner  that  the  results  of  the  pre- 
ceding article  can  be  applied. 

Consider  a  fixed  point 

x  ,y  ,z  and  a  ray  p  ,q  .r  through 
o'o*  0  -o'  o  0     ^^ 


this  point.   Let  (x,y,z)  be  any  point 
on  this  ray.  x,y,  and  z  are  given  by 


(7.1) 


X    =    X 


y  =  y0- 


n 
Jt 


Figure  31 


0         n         0 
For  increasing  Jl  t   the  point   (x,y,z)   travels  on  this  ray  twoards  infinity. 

We  wish  to  determine  the  asymptotic  behavior  of  the  vectors  u",v" 
for  large  values  of  J-  .      We  consider  only  the  vector 


r-  -     BJf  flPil>     e^CW+px+oyH-rz)   dpdq  # 


(7.2)  2^, 

which  is  one  of  the   integrals    (U.15).    We  introduce  the  formulas   (7.1)  into  this 
integral  and  using  the  notation 

(7.3)  r*  =  nTeik(W+pxo  ♦  Wo  +  "0)       , 
it  follows  that 


^u     =   - 


or  letting 

(7.*0 


2tm  SS         r 


-ik  ;::::-(pp  +  qq  +  rr   ) 
n     ^-^o     ^Ho         0 


dpdq 


k«  = 


^^T=  _    |i£l  Jj£  e-i*'(pPo+  qq0  +  rr0) 
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Compsrison  of  this  integral  with  the  integral  in  (7.2)  shows  that 
(7.U)  has  the  same  mathematical  form  as  the  original  diffraction  integral.   Its 
asymptotic  "behavior  for  large  values  of  k'  thus  oust  he  the  same  as  that  of  the 
diffraction  integrals  (U.15)  or  (7.2)  in  case  of  large  values  of  k.   Letting 

(7.5)  <P  =  -    (p?0+  qqQ+  rrQ) 

We  recognize   that  0  has  a  minimum  point  at  p  =  p   ,    q  =  a,   r=r   ;    the  minimum  value 

"being  (p  =  -  n     .     The  results  of  article  6,    formulas    (6.23)  and   (6.6)   in  particular, 
permit  us   to  conclude  that 

X*(po'%)        eik'n2 

(7.6)  IS-    - 


r 


o  ,  /  2 

pp     qq    ^pq 

Lfj£-*~CD    ,  i.e.,  if  k'—i 5-00  .  ihe  quantity,  -  n  ,  plays  the  role  here  that 
-j^*(x,y,z)  plays  as  the  minimum  of  0   in  (7.2).  Now,  at  p  =  p  ,  q  =  q  ,  we  find 
from  (7.5)  that 

P        n 

(7*7)  <*PP  ^qq"  ^pq  =  ^   ' 

o 

Moreover, 

«,  ,    2  ik(W+p  x  *  q  y  *  r  z     -  0n) 

■fa  <  \     ik'n     _       s>7  ,.  vo  0     Ho   0       00* 

A-*   Cp   ,q   )   e  -   n  A(p     q   )      e 


■f,  v      ik>Kx,y,z) 

-     n  Mp   „q   )   e 


o'^o 


for,   hy  formula  (k. 12) 

W(p   ,q    )  +  p  x     +qy     +  r   z       =  ->^"(x.  ,y   ,z    ) 
^o*Mo         ro  0       Ho  0         00  0*    o*    o 

and,    since  V^tx.y, z)   is   the  optical  path  from  the  source  to  (x,y,z), 

">^'(xo'yo'Zo)     ~^n    =V^(x,y,z)   . 
Hence  it   follows  that 
(7.3)  {Lt-T(vQ,%)   eik^    ->    0 

or,    sincei?7f    ->   1,    that 

(7.5)  £    e^^ftp^)      . 

VST  ° 

This,  however,  is  the  condition  (3. 10). 

Ihe  condition  (3.I3)  for  ^he  derivatives  of  u  and  v""can  "be  derived 
in  exactly  the  same  manner. 
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